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Reheating after inflation occurs due to particle production by the oscillating inflaton field. In this 
paper we briefly describe the perturbative approach to reheating, and then concentrate on effects 
beyond the perturbation theory. They are related to the stage of parametric resonance, which we 
called preheating. It may occur in an expanding universe if the initial amplitude of oscillations of 
the inflaton field is large enough. We investigate a simple model of a massive inflaton field cj> coupled 
to another scalar field % with the interaction term g 2 4> 2 X 2 - Parametric resonance in this model is 
very broad. It occurs in a very unusual stochastic manner, which is quite different from parametric 
resonance in the case when the expansion of the universe is neglected. Quantum fields interacting 
with the oscillating inflaton field experience a series of kicks which, because of the rapid expansion 
of the universe, occur with phases uncorrelated to each other. Despite the stochastic nature of the 
process, it leads to exponential growth of fluctuations of the field \. We call this process stochastic 
resonance. We develop the theory of preheating taking into account the expansion of the universe 
and backreaction of produced particles, including the effects of rescattering. This investigation 
extends our previous study of reheating after inflation jlj. We show that the contribution of the 
produced particles to the effective potential V(<j>) is proportional not to cj> 2 , as is usually the case, 
but to \(f>\. The process of preheating can be divided into several distinct stages. In the first stage 
the backreaction of created particles is not important. In the second stage backreaction increases 
the frequency of oscillations of the inflaton field, which makes the process even more efficient than 
before. Then the effects related to scattering of y-particles on the oscillating inflaton field terminate 
the resonance. We calculate the number density of particles n x produced during preheating and 
their quantum fluctuations (x 2 ) with all backreaction effects taken into account. This allows us 
to find the range of masses and coupling constants for which one can have efficient preheating. In 
particular, under certain conditions this process may produce particles with a mass much greater 
than the mass of the inflaton field. 
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I. INTRODUCTION 

According to inflationary theory, (almost) all elemen- 
tary particles populating the universe were created dur- 
ing the process of reheating of the universe after inflation 
. It makes this process extremely important. However, 
for many years the theory of reheating remained the least 
developed part of inflationary theory. Even now, when 
many features of the mechanism of reheating are under- 
stood, the literature on this subject is still full of contra- 
dictory statements. 

The basic idea of reheating after inflation was proposed 
in the first paper on new inflation |||: reheating occurs 
due to particle production by the oscillating scalar field 
<p. In the simplest inflationary models, this field is the 
same inflaton field <f> that drives inflation at the early 
stages of the evolution of the universe. After inflation, 
the scalar field <j> (which we will call inflaton) oscillates 
near the minimum of its effective potential and produces 



elementary particles. These particles interact with each 
other and eventually they come to a state of thermal 
equilibrium at some temperature T. This process com- 
pletes when all (or almost all) the energy of the classical 
scalar field <j> transfers to the thermal energy of elemen- 
tary particles. The temperature of the universe at this 
stage is called the reheating temperature, T r . 

A first attempt at a phcnomenological description of 
this process was made in ref. Q . The authors added var- 
ious friction terms to the equation of motion of the scalar 
field in order to imitate energy transfer from the inflaton 
field to matter. However, it remained unclear what kind 
of terms should be added and whether one should add 
them at the stage of slow rolling of the inflaton field, or 
only at the stage of rapid oscillations of the inflaton field. 

The theory of reheating in application to the new in- 
flation scenario was first developed in refs. and, in 
application to R 2 inflation, in ref. [j7|. It was based on 
perturbation theory, which was quite sufficient for ob- 
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taining the reheating temperature, T r , in many realistic 
models. We will give a detailed description of this theory 
and develop it even further in a forthcoming publication 
|| . However, perturbation theory has certain limitations, 
which have been realized only very recently. In particu- 
lar, the mechanism of decay of the inflaton field to the 
vector fields discussed in || is efficient only at an in- 
termediate stage of reheating in the new inflation model 
considered. The decay of the inflaton field to fermions 
described in || typically is important only at very late 
stages of reheating. In many inflationary models neither 
of these mechanisms gives a correct description of the 
first stages of the process. 

Indeed, recently it was understood Q] that in many 
inflationary models the first stages of reheating occur in 
a regime of a broad parametric resonance. To distin- 
guish this stage from the subsequent stages of slow re- 
heating and thermalization, we called it preheating. The 
energy transfer from the inflaton field to other bose fields 
and particles duringpreheating is extremely efficient. As 
we pointed out in fjj], reheating never completes at the 
stage of parametric resonance; eventually the resonance 
becomes narrow and inefficient, and the final stages of the 
decay of the inflaton field and thermalization of its decay 
products can be described by the elementary theory of re- 
heating J||6j|^]. Thus, the elementary theory of reheating 
proves to be very useful even in the theories where reheat- 
ing begins at the stage of parametric resonance. However, 
it should be applied not to the original coherently oscil- 
lating inflaton field, but to the products of its decay, as 
well as to the part of the inflaton field which survived 
preheating. The short stage of explosively rapid preheat- 
ing in the broad resonance regime may have long-lasting 
effects on the subsequent evolution of the universe. It 
may lead to specific nonthermal phase transitions in the 
early universe (|Jf^l and to topological defect production, 
it may make possible novel mechanisms of baryogenesis 
jn],[l2| , and it may change the final value of the reheating 
temperature T r . 

The theory of parametric resonance in application to 
particle production by oscillating external fields was de- 
veloped more than 20 years ago, see e.g. [f3). The meth- 
ods used in this theory were developed mainly for the 
case of narrow parametric resonance. A first attempt to 
apply this theory to reheating after inflation was made by 
Dolgov and Kirilova [|l4| and by Traschen and Branden- 
berger Jl5| for the narrow resonance regime in the context 
of the new inflation. In it was conjectured that the 
parametric resonance in an expanding universe cannot 
lead to efficient reheating. The authors of Ref. JI1| came 
to an important conclusion that parametric resonance in 
new inflation can be efficient. However, their investiga- 
tion of p arametric resonance was not quite correct, see 
Sec. |V] of this paper. 

In any case, at the moment we do not have any consis- 
tent inflationary models based on the new inflation sce- 
nario. The step towards the general theory of reheating 
in chaotic inflation was rather nontrivial. Indeed, the ef- 



fective potential in new inflation is anomalously flat near 
= 0. As a result of this fine-tuned property of the 
effective potential, the Hubble constant at the end of in- 
flation in this scenario is much smaller than the mass of 
the oscillating scalar field. Therefore the effects related 
to the expansion of the universe are not very destruc- 
tive for the development of the resonance, which may be 
rather efficient even if the resonance is narrow. Narrow 
resonance can be rather efficient in chaotic inflation as 
well, in the context of conformally-invariant theories of 
the type of A^ 4 . In such theories the expansion of the 
universe does not interfere with the development of the 
resonance, and therefore preheating may be efficient even 
if the resonance is rather narrow p|,p^I9| . However, gen- 
erally the effective potential is quadratic with respect to 
<f> near the minimum of the potential, which breaks the 
conformal invariance. As we will show in this paper, for 
the simplest models of inflation, such as the theory of a 
massive inflaton field (j) with quadratic effective poten- 
tial and interaction g 2 4> 2 "x 2 , preheating is efficient only 
if the resonance is extremely broad. The theory of a 
broad parametric resonance in an expanding universe is 
dramatically different from the theory of a narrow reso- 
nance. 

The basic features of the theory of a broad parametric 
resonance were outlined in jl]], where the theory of pre- 
heating was developed in the context of the chaotic in- 
flation scenario, taking into account backreaction of cre- 
ated particles and the expansion of the universe. This 
issue was studied later by many other authors, and a lot 
of very interesting results on parametric resonance and 
particle production have been obtained jf6| - |H| . Of all 
these papers one is especially relevant to our investiga- 
tion. Khlcbnikov and Tkachev |50| performed a detailed 
three-dimensional numerical lattice simulation of broad 
parametric resonance in an expanding universe, taking 
into account the backreaction of produced particles, in- 
cluding, in particular, their rescattering. Their method 
(see also [^J,^8[^9)) is based on solving numerically the 
classical equations for fluctuations of all interacting fields. 
It is presumably the best way to perform computer sim- 
ulations of preheating. 

From the point of view of analytical investigation of 
preheating in the broad resonance regime we should men- 
tion ref. pjl |, where this regime was investigated for the 
case of a non-expanding universe, and some of the re- 
sults of ref. jjj concerning this regime were obtained by 
a different method. However, after our paper Jl| there 
was not much progress in analytical investigation of the 
broad resonance regime in an expanding universe. This is 
not very surprising, because the analytical investigation 
of preheating including backreaction is very difficult; one 
must describe a system of particles far away from equi- 
librium in the regime where effective coupling becomes 
strong because of anomalously large occupation numbers 
of bose particles produced by parametric resonance. But 
the main problem was related to the very unusual na- 
ture of broad parametric resonance in an expanding uni- 
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verse. As we will show in this paper, instead of staying 
in a particular resonance band, each growing mode scans 
many stability/instability bands within a single oscilla- 
tion of the inflaton field, so the usual concept of separate 
resonance bands becomes inadequate. It is a stochastic 
process, during which the number of produced particles 
changes in a chaotic way. On average, the number of 
produced particles grows exponentially, but at some mo- 
ments their number may decrease; a process which would 
be impossible at the classical level. We call this process 
stochastic resonance. The standard methods developed 
for investigation of parametric resonance simply do not 
apply here, so it was necessary to develop a new, more 
general approach. 

The main purpose of the present paper is to develop 
the theory of preheating with an account taken of the 
expansion of the universe and the backreaction of created 
particles, including the effects of their rescattering. We 
will give here a detailed derivation of the results of Rcf. 
0, and describe recent progress in the understanding 
of physical processes which occur soon after the end of 
inflation. 

We will begin our paper with discussion of the evo- 
lution of the scalar fields after inflation neglecting the 
effects of reheating, see Sec. H Sec. Ill contains an in- 
troduction to the elementary theory of reheating |^,||,D . 
We will then develop the theory of particle production 
due to parametric resonance following yj] . First of all, in 
Sec. IV we introduce the theory of reheating due to para- 
metric resonance and discuss the relation between this 
theory and the elementary theory of reheating. Both the- 
ories are very simple, but the transition from one to the 
other is quite nontrivial; it is very difficult to understand 
the theory of parametric resonance using the elementary 
theory of reheating as a starting point, and, conversely, 
perturbation theory is not simply a limiting case of a 
weak parametric resonance. A more detailed discussion 
of all these issues will be contained in our forthcoming 
paper §. 

In Sec. we discuss the difference between the nar- 



Sec. VIII contains a discussion of the backreaction 



row and broad resonance regimes. Sec. VI is devoted 



to a qualitative description of the development of broad 
resonance in an expanding universe. We describe the ef- 
fect of stochastic resonance and illustrate this effect by 
solving the resonance equations numerically, taking into 
account the expansion of the universe. We find that it is 
much easier to perform the investigation in terms of the 
number of created particles, which is an adiabatic invari- 
ant, rather than in terms of wildly oscillating quantities 
such as (x 2 ) which arc studied in many publications on 
preheating. In particular, in some cases (x 2 ) continues 
to grow even after the resonance ceases to exist and the 
number of x particles remains constant. In Sec. VII 



we develop analytic methods for the description of broad 
resonance. These methods are especially appropriate for 
the investigation of stochastic resonance. They are appli- 
cable in those cases where the standard approach based 
on the investigation of Mathicu or Lame equations fails. 



of the x-P ar ti c l es created by parametric resonance on 
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the effective potential of the inflaton field. In Sec 
we describe the process of reheating in the broad res- 
onance regime with an account taken of the change of 
the frequency of oscillations of the inflaton field due to 
its interaction with the x _ P ar ticles produced during pre- 
heating. In Sec. [x| we discuss the process of rescattering 
of X"P ar ticles and the production of 0-particles in this 
process. We also consider some modifications of the pic- 
ture of the second stage of reheating with an account 
taken of rescattering. We calculate the number of par- 
ticles produced during reheating and the amplitude of 
perturbations (x 2 )- In Sec. XI we investigate the possi- 
bility of a copious production of particles with ma ss mu ch 
greater than the inflaton mass. Finally, in Sec. XII we 



give a summary of our results and discuss their possible 
implications. 



II. EVOLUTION OF THE INFLATON FIELD 

During inflation the leading contribution to the energy- 
momentum tensor is given by the inflaton scalar field <fi 
with the Lagrangian 



(i) 



where V{4>) is the effective potential of the scalar field <f>. 
The evolution of the (flat) FRW universe is described by 
the Einstein equation 



H 2 



8tt 
3~Mf I 2 



(2) 



where H = a/a. The Klein-Gordon equation for <p(t) is 



3F< 



(3) 



For sufficiently large initial values of <j> > M p , the "fric- 
tion" term 3H <j> in (||) dominates over (f> and the poten- 
tial term in (^) dominates over the kinetic term. This is 
the inflationary stage, where the universe expands quasi- 
exponentially, a(t) = aoexp(/ dtH(t)). For definiteness, 
we will consider here the simplest models of chaotic infla- 
tion: V((f>) = ^mcj) 2 Q|. In these models inflation occurs 
at <fi > M p . Density perturbations responsible for large- 
scale structure formation in these models are produced 
at <fi ~ 3 — 4 M p . With a decrease of the field <f> be- 
low Mp the "friction" term 3H <f> becomes less and less 
important, and inflation terminates at 4> ~ M p /2. 

When making numerical estimates one should take into 
account that at the last stages of inflation the friction 
term is still non-negligible, and therefore during the first 
oscillation the amplitude of the field rapidly drops down. 
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FIG. 1. Oscillations of the field <j> after inflation in the 

2,2 

theory m ^ . The value of the scalar field here and in all 
other figures in this paper is measured in units of M p , time is 
measured in units of m _1 . 



For the quadratic potential V(tf>) 



the ampli- 



tude after the first oscillation becomes only 0.04M p , i.e. 
it drops by a factor of ten during the first oscillation, 
see Fig. 1. Later on the solution for the scalar field 4> 
asymptotically approaches the regime 



4>(t) = $(t) • sin mi , 



*(*) = 



v 



3nmt 2vrV37riV 



(4) 



Here $(£) is the amplitude of oscillations, N is the num- 
ber of oscillations since the end of inflation. For simple 
estimates which we will make later one may use 



*(*) 



M p 
3mt 



M p 
20N' 



(5) 



The scale factor averaged over several oscillations grows 
as a(t) w ao(-^) 2 ^ 3 . Oscillations of <p in this theory 
are sinusoidal, with the decreasing amplitude $(t) = 

^lf~(~a(t)) 3 ^ 2 ■ ^ ne ener §y density of the field <f> decreases 
in the same way as the density of nonrelativistic parti- 



cles of mass 



3 . Hence the 



coherent oscillations of the homogeneous scalar field cor- 
respond to the matter dominated effective equation of 
state with vanishing pressure. 

Reheating occurs when the amplitude of oscillations of 
the inflaton field <fi decreases much faster than in (Eh , and 
its energy density is transferred to the energy density of 
other particles and fields. 



III. OSCILLATIONS AND DECAY OF THE 
SCALAR FIELD 

In the present section, we will discuss the elementary 
theory of reheating developed in [|||]; see also j|. A 
more detailed discussion of this theory will be contained 
in M. We will consider a basic model describing the 



inflaton scalar field i 
and a spinor field i/>: 



interacting with a scalar field x 



L 



V(4>) 



2 /2,,2 



(0) x 2 + ^Rx 2 
- h<ijnp<t> ■ ( 6 ) 



Here g, h and £ are small coupling constants, R is the 
space-time curvature, and V(<p) is the effective potential 
of the field <fi. We will suppose here, for generality, that 
the effective potential has a minimum at <f> = a, and 
near the minimum it is quadratic with respect to the 



field 0: V((t>) 



2 {4> — a) 2 . Here m 2 is the effective 
mass squared of the field <p. After the shift <fi ~ a ~~ * 4>i 
the effective potential acquires the familiar form ^m 2 (j) 2 , 
and the Lagrangian acquires an interaction term which is 
linear with respect to the field (j>: AL = —g 2 a4>x 2 - This 
term vanishes in the case without spontaneous symmetry 
breaking, where a = 0. The masses of the x-particles 

and ip after the shift become m x — ^m 2 (0) + g 2 a 2 and 

= m^,(0) + ha. In this section we will consider the 
case m ^> m x , m^. We will assume that after inflation 
H -C m. This condition is always satisfied during the 
last, most important stages of reheating. 

We will study now the oscillation of the scalar field 
near the minimum of its effective potential. The energy 
density of the oscillating field (after the shift <f> — a — > <f>) 
is = j4> 2 + ^m 2 cf> 2 . If we ignore for a moment the 
effects associated with particle creation, the field <fi af- 
ter inflation oscillates near the point <fi = with the 
frequency fco = m. The amplitude of oscillation de- 
creases as a~ 2 due to the expansion of the universe, 
and the energy of the field <j> decreases in the same way 
as the density of nonrelativistic particles of mass m: 
Pcf, = \4> 2 + \m 2 cf) 2 ~ a~ 3 . A homogeneous scalar field 
oscillating with frequency m can be considered as a co- 
herent wave of (^-particles with zero momenta and with 
particle density — p^/m. In other words, oscilla- 
tors of the same frequency m, oscillating coherently with 
the same phase, can be described as a single homoge- 
neous wave <j)(t). Note that if we consider time intervals 
larger than the typical oscillation time to -1 , the energy 
density of the oscillating field, and the number density of 
the particles n$ will be related to its amplitude $ in a 
simple way: 



(7) 



(8) 



Now we will consider effects related to the expansion of 
the universe and to particle production. For a homoge- 
neous scalar field in a universe with a Hubble constant H, 
the equation of motion with non-gravitational quantum 
corrections is 



(f> + 3H(t) <j)+(m 2 + U(u>)) = 0. 



(9) 
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Here H(uj) is the flat space polarization operator for the 
field <f> with four-momentum fcj = (w, 0, 0, 0), ui = to. 

The real part of II (a; ) gives only a small correction 
to to 2 , but when lo > min(2m x , 2m^) , the polariza- 
tion operator TL(u>) acquires an imaginary part lmll(uj). 



We will assume that 



» H l 



> Imn. The 



first condition is automatically satisfied after the end 
of inflation; the second is usually also true. We have 
= $ a" 3/2 0) = $ exp(-§ f dtH(t)). Neglecting 
for simplicity the time-dependence of H and Im II due to 
the expansion of the universe, we obtain a solution of (j^) 
that describes damped oscillations of the field near the 
point = 0: 



<fi = 3>(t) exp(iTOi) 
4>q exp(iTOi) • exp 



--13H 



ImlT(TO) 



From unitarity it follows that |B2[ 



Imll = mT 



(10) 



(11) 



where T = T(4> — > xx) + F(0 — > is the total decay 
rate of (^-particles. (In a more general case one should 
calculate not only the imaginary part of the polarization 
operator, but the imaginary part of the effective action 
0.) Thus, Eq. ( |l(i| ) implies that the amplitude of os- 
cillations of the field <f> decreases as exp 
due to particle production which occurs during the decay 
of the inflaton field. 

Note that under the condition m 3> H, the polariza- 
tion operator II and the decay rates T do not depend 
on the curvature of the universe (and thus on time) and 
coincide with their flat-space limits. In particular, the 
probability of decay of a 0-particle into a pair of scalar 
X-particles or spinor ^-particles for m 3> m x , is given 
by the following expressions : 



r(</> -> xx) 



87TTO 



Sir 



(12) 



For a phenomenological description of the damping of 
oscillations of the scalar field <f> ( [Toj ) one may add an 
extra friction term T<f> to the classical equation of motion 
of the field (p, instead of adding the term proportional to 
the imaginary part of the polarization operator, 



<f> + 3H(t) <f> + T<f> + m 2 (j) = 



(13) 



This phenomenological equation together with relation 
( |ll| ) for T reproduces the damped oscillator solution jic| ) 
of Eq. (g). The idea that one can describe effects of 
reheating by adding friction terms to the equation of 
motion goes back to one of the first papers on reheat- 
ing jj]. At first the physical origin of such terms as well 
as their value remained obscure. Some authors added 
various auxiliary friction terms to the equations of the 
inflaton field in order to slow down its motion and make 
inflation longer, see e.g. 0J33[. From the derivation of 



expression ([11]) for Y it follows, however, that the simple 
phenomenological equation (|13|) is valid only at the stage 
of rapid oscillations of the field <j> near the minimum of 
V(<p). This equation cannot be used to investigate the 
stage of slow rolling of the field <j> during inflation. 

According to (J10|), the field amplitude <!>(£) obeys the 
equation 



1 d 



: (a 3 $ 2 ) = -r$ 2 



(14) 



If one multiplies it by to, one obtains the following equa- 
tion for the number density @ of the coherently oscil- 
lating 0-particles: 



d_ 
~dt 



(a 3 ^) = -r ■ a 3 



(15) 



This equation has a simple interpretation. It shows that 
the total comoving number density of particles ~ a^n^ 
exponentially decreases with the decay rate T. Similarly, 
one obtains the following equation for the total energy of 
the oscillating field <fi: 



= -T ■ a A 



(16) 



The decay products of the scalar field (j> are ultrarela- 
tivistic (for to ^> m x ,m^), and their energy density de- 
creases due to the expansion of the universe much faster 
than the energy of the oscillating field (p. Therefore, re- 
heating in our model ends only when the Hubble con- 
stant H ~ H becomes smaller than T, because otherwise 
the main portion of energy remains stored in the field <fr. 
Therefore the age of the universe when reheating com- 
pletes is given by t r ~ fr _1 - At that stage the main 
part of the matter in the universe becomes ultrarelativis- 
tic. The age of the universe with the energy density p is 
t = V~Mp/ \f&Tip Q. This, together with the condition 
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gives the energy density at time t r 



p(t r 



3r 2 Af p 2 
8tt 



(17) 



If thermodynamic equilibrium sets in quickly after the 
decay of the inflaton field, then the matter acquires a 
temperature T r , which is defined by the equation 



p(t r ) 



3r 2 Af 2 
8tt 



7r 2 iV(T r ) T4 



30 



(18) 



Here N(T) is the number of relativistic degrees of free- 
dom at the temperature T; one should take 1 for each 
scalar, two for each massless vector particle, etc. ||. In 
realistic models one may expect N(T r ) ~ 10 2 — 10 3 , which 
gives the following estimate of the reheating temperature: 



(19) 



Note that T r does not depend on the initial value of the 
field 4>; it is completely determined by the parameters of 
the underlying elementary particle theory. 
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Here we should make an important comment. In the 
absence of fermions, the only contribution to the decay 

4 2 

rate would be — ► xx) = fi^r- Note that this term 
disappears in the theories without spontaneous symme- 
try breaking, where a — 0. This does not necessarily 
mean that there is no reheating at all in such theories. 
Indeed, decay is possible not only in the presence of a con- 
stant field a but in the presence of a large oscillating field 
4>{t) as well. In what follows we will study parametric res- 
onance and reheating in models with a — 0, or a <C $, 
where $ is the amplitude of the oscillations. However, 
when reheating proceeds and $ becomes small one may 
expect perturbation theory to work well. To get an esti- 
mate for the decay rate at a — let us simply write $ 
instead of a in Eq. (|l|): Y{<j)(f) -> X x) ~ f^- The prob- 
lem with this term is that $ 2 decreases as t~ 2 in the ex- 
panding universe, whereas the Hubble constant decreases 
only as f -1 . Therefore the decay rate never catches up 
with the expansion of the universe, and reheating never 
completes. Reheating can be complete only if T decreases 
more slowly than t . Typically this requires either spon- 
taneous symmetry breaking (er ^ 0) or coupling of the 
inflaton field to fermions with < m/2. If both of 
these conditions are violated, the inflaton field <j) never 
decays completely. Such fields may be responsible for the 
dark matter of the universe, but it requires certain fine- 
tuning of the parameters. Normally, an incomplete decay 
of the inflaton field implies that the universe at the age 
of 10 billion years is cold, empty and unsuitable for life. 
We should emphasize that this may happen even if the 
coupling constant g 2 is very large. Thus the requirement 
that reheating is complete imposes important constraints 
on the structure of the theory. 

The elementary theory of reheating described above is 
simple and intuitively appealing. It proves to be very 
successful in describing reheating after inflation in many 
realistic inflationary models. That is why we are going 
to develop this theory even further in However, in 
some cases where the amplitude of the oscillating field is 
sufficiently large, reheating occurs in a different way, in 
the regime of parametric resonance. 



IV. PARAMETRIC RESONANCE AND LIMITS 
OF APPLICABILITY OF PERTURBATION 
THEORY 

A. Perturbation theory versus narrow resonance 

In the investigation performed above we made a natu- 
ral assumption that the decay probability T of the scalar 
field 4> can be calculated by ordinary methods of quan- 
tum field theory describing the decay (f> — > XX- However, 
if many x-particles have already been produced, rife > 1, 
then the probability of decay becomes greatly enhanced 
due to effects related to Bose-statistics. This may lead 
to explosive particle production. 



For simplicity, we consider here the interaction be- 
tween the classical inflaton field <f> and the quantum 
scalar field x with the Lagrangian (^|). The Heisenberg 
representation of the quantum scalar field x is 



f 



(2tt) 3 / 2 



d A k 



(a k Xk(t)e-^ + a+xl(t), 



(20) 



where and are annihilation and creation operators. 
For a flat Friedmann background with scale factor a(t) 
the temporal part of the eigenfunction with comoving 
momentum k obeys the following equation: 



Xk + 3-x fc + -J + r<(0) - £R + 5 V )Xk = . (21) 



(The physical momentum p = jjfer coincides with k for 
Minkowski space, where a = 1.) Eq. ( pl| ) describes an 
oscillator with a variable frequency u due to the time- 
dependence of a(t) and the background field <p(i). Un- 
til the last section of this paper we will suppose that 
the effective mass of the field x vanishes for 4> = 0: 
m x (0) = 0. 



In Sec. XI we will investigate the opposite 



case, m x (0) ^> m. 

As in the previous section, consider the simplest po- 
tential V{4>) ~ ±m 2 (4> - a) 2 (to mimic the situation 
with spontaneous symmetry breaking) and make the shift 
<f> — g —> cj), after which the effective potential becomes 
im 2 </) 2 , and the interaction term — \g 2 (j) 2 x 2 transforms 
to —\g 2 4> 2 x 2 ^ q 2 Q(t>X 2 — \g 2 o 2 x 2 ■ We shall analyze the 
general equation ( |2l|) in different regimes. 

Suppose first that the amplitude of oscillations <fi is 
much smaller than <j, and neglect for a moment the ex- 
pansion of the universe, taking a = 1 in Eq. (pl|). Then 
one can write the equation for modes (quantum fluctu- 
ations) of the field x with physical momentum k in the 
following form: 



Xk 



(k 2 



2 2 

go 



2g cr$ sin mi) Xk = 



(22) 



where k = vk 2 , and <I> stands for the amplitude of oscil- 
lations of the inflaton field. This equation describes an 
oscillator with a periodically changing frequency lo 2 (t) — 
k 2 + g 2 a 2 + 2g 2 cr$ sinmt. This periodicity may lead to 
parametric resonance for modes with certain values of k. 
The simplest way to describe this important effect is to 
make a change of variables mt = 2z — tt/2, which reduces 
Eq. (G2) to the well-known Mathieu equation |34] : 



Xk + (A k -2qcos2z)xk = 



(23) 



Here Ai 



i 2 i 2 2 



= £2!£* z = irU d ■ 
rti 2 ' 2 5 1 

denotes differentiation with respect to z. The properties 
of the solutions of the Mathieu equation are well rep- 
resented by its stability/instability chart which can be 
found, e.g., in B4]. An important feature of solutions 



G 



of Eq. ( J23| ) is the existence of an exponential instabil- 
ity Xk oc exp(/i^"^z) within the set of resonance bands of 
frequencies Ak^ labeled by an integer index n. This in- 
stability corresponds to exponential growth of occupation 

numbers of quantum fluctuations cx exp(2/zjj™^) 

that may be interpreted as particle production. In a state 
with a large number of Bose particles the estimates for F 
obtained in the previous subsection do not apply, and one 
should use much more elaborate methods of investigation 
based on the theory of parametric resonance. 

In the case under consideration, g$ <C go -C m, the 
theory of parametric resonance is well known p5[ . In- 
deed, in this case one has q <^ 1, and the resonance oc- 
curs only in some narrow bands near A k ~Z 2 , I = 1,2, .... 
Each band in momentum space has width of order Afc ~ 
q l , so for q < 1 the widest and most important instability 
band is the first one, Ak~l±q = l± ig_gg _ 

The factor [i k which describes the rate of exponential 
growth for the first instability band for m 2 3> g 2 o 2 is 
given by @ 




Thus resonance occurs for k = ?r(l ± I). The index 
vanishes at the edges of the resonance band and takes 
its maximal value u k = % = 2 ^- at k = S. The 

' K 2 m z 2 

corresponding modes Xk grow at a maximal rate exp(^), 

which in our case is given by exp( 2 ^) = exp( g < ^ t ). 

The growth of the modes Xfc leads to the growth of 
the occupation numbers of the created particles n k (t). 
Indeed, the number density n k of particles with mo- 
mentum k can be evaluated as the energy of that mode 
^IXfel 2 + \u 2 \xk\ 2 divided by the energy uj k of each par- 
ticle: 

(\Xk\ 2 , 2 \ 1 fOK\ 

nk = Y{^f + lxkl J' 2- (25) 

When the modes Xfc grow as exp(^), the number of 
X-particles grows as exp(gz), which in our case is equal 
toexp(2f*)=exp(2^*l). 

The fact that the resonance occurs near k = ^ has a 
simple interpretation: In the limit go <C m the effective 
mass of the x-particles is much smaller than m. There- 
fore one decaying (^-particle creates two x-particles with 
momentum k ~ m/2. This picture is very similar to 
the process of decay (f> — > XX discussed in the previous 
section, but the results are absolutely different. Indeed, 
in perturbation theory the amount of produced parti- 
cles did not depend on the number of particles produced 

earlier, and the rate of production for our model was 

4 2 

given by T((j) — > xx) — smn ' Thus the decay rate T 
was suppressed by the factor g 4 , which made the decay 
very slow in the weak coupling limit. By contrast, in the 
regime of parametric resonance the rate of production of 



X-particles is proportional to the amount of particles pro- 
duced earlier (which is why we have exponential growth) , 
and the rate of the process is given by an absolutely dif- 
ferent expression /ifem ~ 9 ^ , which is greater than T 
for $ > (-0. 

Thus, before going any further we should understand 
how these two processes are related to each other, and 
why we did not find the effect of parametric resonance 
in the investigation performed in the previous section. Is 
the perturbation theory discussed there a limiting case 
of the narrow resonance regime or is it something else? 

The reason we missed the effect of parametric reso- 
nance is rather delicate. In our calculations of the imag- 
inary part of the polarization operator we assumed that 
the x-particles produced by the oscillating scalar field <f> 
are normal particles on the mass shell, fc£ = m 2 , . This is 
what one would get solving Eq. (p2|) in any finite order of 
perturbation theory with respect to the interaction term 
2g 2 xo^ sinmt. However, if one solves the equation for 
the fluctuations of the field x ( ^2|) exactly, one finds expo- 
nentially growing modes Xfe- This creates a new channel 
of decay of the scalar field 4>. 

Note that exponentially growing modes occupy a very 
small portion of momentum space in the narrow reso- 
nance limit. This means that the fluctuations of the field 
X for almost all k are normal fluctuations with k 2 , = m 2 ,. 
In this case our calculation of the imaginary part of the 
polarization operator does apply. If the resulting value 
of r appears to be smaller than ~ qm, then the 

perturbative decay of the scalar field may coexist with 
the parametric resonance. One may consider several dif- 
ferent possibilities. In the beginning the scalar field <p 

2 

can be expected to oscillate with amplitude <& > §^cr. In 
this regime parametric resonance leads to the exponen- 
tial growth of modes Xk, as we discussed above. However, 

gradually the field <j> loses its energy, and its amplitude 

2 

$ becomes smaller than f^cr. In this regime the ampli- 
tude of the field $ decays exponentially within a time T _1 
which is smaller than the typical time necessary for para- 
metric resonance to occur. One may say that the pertur- 
bative decay makes the energy eigenstate (the mass) of 
the field <f> "wide," with width T, and when this width 
exceeds the width of the resonance band ~ qm/2, the res- 
onance terminates. Starting from this moment perturba- 
tion theory takes over, and the description of reheating 
should be given along the line of the elementary theory 
described in the previous section. 

Thus, the standard effect of scalar field decay described 
by the elementary theory of reheating and preheat- 

ing due to parametric resonance are two different effects. 
In an expanding universe there exist other reasons for 
evolving from parametric resonance to perturbative de- 
cay. 

First of all, during the expansion of the universe the 
field (/> decreases not only because of its decay, but be- 
cause of the "friction term" 3H <fi in the equation of mo- 
tion for the field <f>. Thus one should compare qm with the 
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effective decay rate 3iJ + T: Parametric resonance occurs 
only for qm > 3H + T. Note that for T > H perturba- 
tive decay leads to reheating even neglecting parametric 
resonance. Therefore to check whether parametric res- 
onance appears at the time when perturbative decay is 
inefficient, i.e. in the case r < H , it is enough to consider 
the condition qm > 3H. 

Another important mechanism which can prevent 
parametric resonance from being efficient is the redshift 
of momenta k away from the resonance band. The to- 
tal width of the first band is given by qm; the width 
of the part where the resonance is efficient is somewhat 
smaller; one can estimate it as 222. The time At during 
which a given mode remains within this band depends on 
the equation of state of matter, and typically can be esti- 
mated by qH~ x . During this time the number of particles 
in growing modes increases as exp^2_2Ii^ , This leads to 

efficient decay of inflatons only if q 2 m > H . In the nar- 
row resonance limit q <C 1 this is a stronger condition 
than the condition qm > 3H. 

In general, it is possible that exponential growth dur- 
ing the time At is small, but At <C H^ 1 and therefore 
resonance still plays some role in reheating. However, 
this is a rather exceptional situation. Therefore typically 
the set of conditions for the resonance to be efficient can 
be formulated as follows: 



qm > r , 
q 2 m>H 



(26) 
(27) 



In the model considered above these conditions yield: 
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n\JmH 
Ag 2 a 



(28) 
(29) 



Thus parametric resonance can be efficient at a suffi- 
ciently large but reheating never ends in the regime 
of parametric resonance. As soon as the amplitude of 
oscillations becomes sufficiently small, parametric reso- 
nance terminates, and reheating can be described by the 
elementary theory developed in ||,||,|| ■ Typically the re- 
heating temperature is determined by these last stages 
of this process. Therefore one should not calculate the 
reheating temperature simply by finding the endpoint of 
the stage of parametric resonance, as many authors do. 
The role of the stage of preheating is to prepare a differ- 
ent setup for the last stage of reheating. It changes the 
reheating temperature, and it may lead to interesting ef- 
fects such as nonthermal symmetry restoration and new 
mechanisms of baryogenesis. However, reheating never 
ends in the regime of parametric resonance; it does not 
make much sense to calculate the reheating temperature 
at the end of the stage of preheating. 



The expansion of the universe and the inflaton decay 
are not the only mechanisms which could prevent the de- 
velopment of resonance. As we will show, backreaction of 
created particles may change the parameters A% and q. 
There will be no resonance if the x-particles decay with 
decay rate T > /i^m, or if within the time ~ (/x^m) -1 
they are taken away from the resonance band because of 
their interactions. Also, there is no explosive reheating 
if the decay products include fermions since the fermion 
occupation numbers cannot be large because of the Pauli 
principle. This happens, for example, in many inflation- 
ary models based on supergravity where inflaton decay 
is often accompanied by gravitino production . 

If reheating never ends in a state of narrow parametric 
resonance, one may wonder whether reheating may begin 
in a state of narrow resonance. As we are going to show, 
in most cases inflation begins in a state of broad para- 
metric resonance; the resonance typically ceases to exist 
as soon as it becomes narrow. But before analyzing this 
issue, we will take one last look at the model which we 
studied above. 



B. Processes at <j> ~ a 

In our investigation of the simple model with sponta- 
neous symmetry breaking (ct ^ 0) we assumed that the 
amplitude of oscillations of the scalar field <j> is very small, 
$ -C o. Therefore we retained only the quadratic part 
of the effective potential, V(cp) ~ (<j) — ct) 2 . However, in 
realistic models of spontaneous symmetry breaking this 
condition is satisfied only at the end of parametric reso- 
nance. Indeed, let us consider a theory with spontaneous 
symmetry breaking with the usual potential j((j) 2 — a 2 ) 2 . 
Then after spontaneous symmetry breaking and the cor- 
responding shift cf> — a — > (j> the theory at -C a can 
be represented as a theory of a massive scalar field with 
a mass m 2 



2A<7 2 interacting with the field \ which 
acquires mass m 2 = g 2 u 2 . In this respect, it coincides 
with the toy model studied in the previous subsection. 
However, there are some important differences. 

First of all, the process <f> — » XX is possible only if 
m > 2m x . This was one of the conditions which we 
used in our investigation: we assumed that m 3> m x , 
i.e. X ^> g 2 . However, in this case the interaction j(j) 4 
which we did not take into account so far may become 
more important than the interaction \4> 2 x 2 which we 
considered. As a result, the production of ^-particles 
may be more efficient than the production of x-particles. 

In order to investigate this possibility let us study for 
a moment a model with the effective potential j(<j) 2 — 
a 2 ) 2 in the limit X ^S> g 2 , i.e. neglecting the interaction 
%f4> 2 X 2 - We will assume here that in the beginning the 
field <j) was & t the top of the effective potential. At that 
time its effective mass squared was negative, m 2 (0) = 
—Act 2 . This fact alone, independent of any parametric 
resonance, leads to the production of particles of the field 
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cf>. The main point here is that all modes with k < VAc 
grow exponentially, which breaks the homogeneity of the 
oscillating scalar field. This is an interesting effect, which 
has some nontrivial features, especially if one takes the 
expansion of the universe into account. We will return 
to its discussion elsewhere. However, this effect does not 
last long because away from the maximum of the effective 
potential its curvature becomes positive. 

When the amplitude of the oscillations of the field <f> 
near <j) = cr becomes smaller than a, the field begins os- 
cillating near its minimum with a frequency m « \[T\o. 
The parametric resonance with </>-particle production in 
this regime can be qualitatively understood if the equa- 
tion for the fluctuations S4>k is approximately represented 
as a Mathieu equation. The modes 4>k grow in essentially 
the same way as the modes in the second instability band 
of the Mathieu equation with Au = + 4, q — 6^. For 
q JSJ 1, we are in the broad resonance regime, and there 
is a significant production of ^-particles. However, for 
q <C 1, i.e. for $ <C c/6, the parametric resonance in the 
second band becomes very inefficient. (One can obtain 
the same result by a more accurate investigation of para- 
metric resonance in this situation in terms of the Lame 
equation, but this is not our purpose here.) 

Thus, we are coming to the following picture of para- 
metric resonance in this model. In the beginning of the 
rolling of the field <j) down to the minimum of the ef- 
fective potential, the leading source of particle produc- 
tion is associated with the tachyonic mass of the field 
4>. Soon after that, the leading mechanism is the decay 
of a coherently oscillating field cf> into ^-particles. This 
mechanism remains dominant until the amplitude of the 
field $ becomes much smaller than a, after which the 
decay <fi — > XX studied in the previous section becomes 
more important. (This process becomes somewhat more 
complicated if the backreaction of the produced particles 
it taken into account.) Finally, when the amplitude of 

2 

the oscillations $ becomes smaller than Jrrcr, or when 
it becomes smaller than " ,v '" 
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whichever comes first, 



the parametric resonance ceases to exist, and the decay 
4> ~ * XX is described by the elementary theory of reheat- 
ing based on perturbation theory. 

We should note that the x-particle production in this 
model for A ^ g 2 was first studied in jl5|]. However, as 
we just mentioned, at <& ~ a this process is subdominant 
as compared to the ^-particle production, which was not 
studied in [fl5| . The process of x-particle production is 
more efficient than </>-particle production only for $ <C cr. 
In this regime our results differ from those obtained in 
[ p5| by the factor * in the exponent. This difference is 
very significant because it leads to a much less efficient 
reheating, which shuts down as soon as becomes suffi- 
ciently small. 

The models studied in the last two subsections can be 
considered as a good laboratory where one can study dif- 
ferent features of parametric resonance. However, in our 
investigation so far we did not discuss the question of 



initial conditions for resonance in these models. Indeed, 
after 15 years of investigation we still have not found any 
simple mechanism which will put the inflaton field on the 
top of the potential at (f> = in the new inflation scenario. 
Also, the shape of the potential required for new infla- 
tion (extremely flat near the origin) is rather artificial. 
As soon as we consider generic initial conditions for the 
scalar field 4> in more realistic inflationary models, such 
as chaotic inflation in the theory with a simple potential 

2.2 

m ® , the theory of parametric resonance becomes dif- 
ferent in many respects from the simple theory described 
above. 



V. BROAD RESONANCE VERSUS NARROW 
RESONANCE IN MINKOWSKI SPACE 

In the chaotic inflation scenario one does not impose 
any a priori conditions on the initial value of the scalar 
field. In many models of chaotic inflation the initial am- 
plitude of oscillations of the field <fi can be as large as 
M p , i.e. much greater than any other dimensional pa- 
rameters such as cr. Therefore in the remaining part of 
the paper we will concentrate on the simplest chaotic 
inflation model without symmetry breaking with the ef- 
fective potential V(cb) = ■ r ^-^ 2 , and the interaction term 
— ^g 2 cf> 2 X 2 - In this case instead of Eq. ( p2[ ) one has 

Xk + {k 2 + g 2 <S> 2 sin 2 (rrrf)) X k = . (30) 

This equation describes an oscillator with a periodically 
changing frequency uj 2 (t) = k 2 + g 2 <& 2 sin 2 mt. One can 
write it as a Mathieu equation (Eq. (p3|)) with Ak — 

For <?$ < m we have a narrow resonance with g« 1. 
In this regime the resonance is more pronounced in the 
first resonance band, for modes with k 2 ~ m 2 (l — 2q±q). 
The modes Xk with momenta corresponding to the cen- 

m grow as e qz ^ 2 , which in 

^ and the number 



'( 9 8m 



tcr of the resonance at k 
our case equals e Mfem * ~ exp( 

of x-particles grows as e 2fJ,kZ ~ e qz ~ exp ^ 9 ® m 1 j . This 

process can be interpreted as a resonance with decay of 
particles with mass m to two x _ P ar ticles with mo- 
rn. We show the results of the numerical so- 



two 

menta k 

lution of Eq. (30) for the fastest growing mode Xk in the 
narrow resonance regime in Fig. 2. Typically, the rate of 
development of the parametric resonance does not differ 
much from the rate of the growth of the leading mode 
Xk> see a discussion of this issue in the next section. 
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FIG. 2. Narrow parametric resonance for the field \ in the 

2 ,2 

theory „ 9 in Minkowski space for g ~ 0.1. Time is shown 
in units of m/2n, which is equal to the number of oscillations 
of the inflaton field cf>. For each oscillation of the field <f>(t) the 
growing modes of the field x oscillate one time. The upper 
figure shows the growth of the mode Xfe f° r the momentum 
k corresponding to the maximal speed of growth. The lower 
figure shows the logarithm of the occupation number of par- 
ticles rik in this mode, see Eq. (p^). As we see, the number 
of particles grows exponentially, and Inn* in the narrow res- 
onance regime looks like a straight line with a constant slope. 
This slope divided by 4tt gives the value of the parameter fik- 
In this particular case fik ~ 0.05, exactly as it should be in 
accordance with the relation fik ~ <?/2 for this model. 

On the other hand, for oscillations with a large ampli- 
tude <E> the parameter q = ^J^ 2 can be very large. In this 
regime the resonance occurs for a broad range of values 
of k, the parameter fj, k can be rather large, and reheat- 
ing becomes extremely efficient. The resonance occurs 
for modes with = A — 2q, i.e. above the line A = 2q 
on the stability /instability chart for the Mathieu equa- 
tion . The standard methods of investigation of narrow 
parametric resonance do not work here. The difference 
between these two regimes can be easily grasped by com- 
paring solutions of Eq. (53) for small and for large q, see 
Figs. 2 and 3. 



FIG. 3. Broad parametric resonance for the field \ m 
Minkowski space for q ~ 2 x 10 2 in the theory — ^ . For each 
oscillation of the field <f>(i) the field \k oscillates many times. 
Each peak in the amplitude of the oscillations of the field \ 
corresponds to a place where <f>(t) = 0. At this time the oc- 
cupation number rik is not well defined, but soon after that 
time it stabilizes at a new, higher level, and remains constant 
until the next jump. A comparison of the two parts of this 
figure demonstrates the importance of using proper variables 
for the description of preheating. Both Xk and the integrated 
dispersion {x 2 ) behave erratically in the process of parametric 
resonance. Meanwhile is an adiabatic invariant. Therefore, 
the behavior of nt is relatively simple and predictable every- 
where except the short intervals of time when <j>{t) is very 
small and the particle production occurs. In our particular 
case, the average rate of growth of nk is close to the maximal 
possible rate for our model, \ik ~ 0.3. 

The time evolution is shown in units to/27t, which cor- 
responds to the number of oscillations N of the inflaton 
field 4>. The oscillating field cf>(t) ~ $ sin mi is zero at in- 
teger and half-integer values of the variable mt/2n. This 
allows us to identify particle production with time inter- 
vals when <f)(t) is very small. 

During each oscillation of the inflaton field cf>, the 
field x oscillates many times. Indeed, the effective mass 
m x (t) = g4>(t) is much greater than the inflaton mass m 
for the main part of the period of oscillation of the field 
4> in the broad resonance regime with q 1 / 2 = 3> 1. 
As a result, the typical frequency of oscillation uj(t) = 
yjk 2 + g 2 (f> 2 (t) of the field x is much higher than that of 
the field <f>. Within one period of oscillation of the field <f> 
the field x makes 0(q 1 ^ 2 ) oscillations. That is why dur- 
ing the most of this interval it is possible to talk about 
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an adiabatically changing effective mass m x (t). There- 
fore, in the broad resonance regime the amplitude of \k 

is minimal at the points where the frequency is maxi- 
_ 1 

mal, |xfc| oc bj(t) 2 5 i.e. at <fi(t) = and it increases 
substantially near the points at which <f>(t) = 0. 

For very small <f>(t) the change in the frequency of os- 
cillations w(t) ceases to be adiabatic. The standard con- 
dition necessary for particle production is the absence of 
adiabaticity in the change of w(t): 



~dt 



>cj 2 



(31) 



One should note that for a narrow resonance this condi- 
tion is not necessary, because even a small variation of 
u>(t) may be exponentially accumulated in the course of 
time. However, for a broad resonance one should expect 
a considerable effect du ring each oscillation, which im- 
plies that the condition (|3l| ) should be satisfied. To find 
the time interval At* and the typical momenta fc* when 
and where it may happen let us remember that for small 
(j) one has 4> ~ m$. Therefore our condition (|T]) implies 
that 



fc 2 < ( 5 2 0m$) 2 / 3 



2 ±2 

9 <f> 



(32) 



Let us consider those momenta fc 2 which satisfy condition 
( |32| ) as a function of <j)(t). This condition becomes satis- 
fied for small fc when the field <p(t) becomes smaller than 

'2J£, The maximal range of momenta for which particle 



Each time the field <j> approaches the point <j> = 0, it 
spends time 



At, 



2</>* 



1 



K 1 



(35) 



< 



During that time fc* 



production occurs corresponds to 4>{t) 



where 



in the domain 
g4>*, so that to ~ fc*. This estimate of At* tells us that 
particle production in the broad resonance regime occurs 
within a time of order of the period of one oscillation of 
the field \, At* ~ oj^ 1 , in agreement with the uncertainty 
principle. One can easily identify these short intervals in 
Fig. | 

In the semiclassical regime when the frequency Uk(t) is 
changing adiabatically, rik is a constant which coincides 
with an adiabatic invariant. To appreciate the useful- 
ness of the introduction of the adiabatic invariant rife, 
one should compare the evolution of the modes Xk with 
the evolution of the occupation numbers corresponding 
to each of these modes shown in Figs. || and [| As 
we see, in the narrow resonance regime Xfc vigorously 
oscillates, whereas lnn/c grows like a straight line. In 
the broad resonance regime the field amplification occurs 
near the points 4>(t) — where the process is not adia- 
batic. The occupation number rife, being an adiabatic 
invariant, changes only during these short time intervals, 
when the number of particles is not well defined. 

Analytical investigation of the broad resonance regime 
in the context of the theory of reheating was first reported 
in H , see also . Now we are going to perform a much 
more detailed investigation of this regime. 



1 / m<I> 1 



1/4 



(33) 



VI. STOCHASTIC RESONANCE IN AN 
EXPANDING UNIVERSE 



The maximal value of momentum for particles produced 
at that epoch can be estimated by fc max = y^§^- In the 

main part of the interval \<f>\ < 20* the range of momen- 
tum remains smaller but the same order of magnitude as 
fcmax- Thus one may estimate a typical value of momen- 
tum of particles produced at that stage as fc*/2, where 



fc* = y/gm® = V2 



mq 



1/4 



(34) 



This simple estimate practically coincides with the re- 
sult of a more detailed and rigorous investigation which 
will be performed in Sec. VII. This is a very important 
result 0, which we are going to use throughout the pa- 
per This result implies, in particular, that in the broad 
resonance regime m <C fc* *C g$. 



*In this paper we will use both physical momenta and co- 
moving momenta. Our definition of fc* refers to physical 
momentum. 



To understand why the broad resonance regime is so 
important for the theory of reheating in an expanding 
universe, let us remember that resonance in an expanding 
universe appears only for q 2 m > H , which in our case 
reads 



g<I> > 2m 



H 



1/4 



(36) 



In the simplest inflationary models including the model 
which we consider now the value of the Hubble constant 
at the end of inflation is of the same order as the inflaton 
mass m, but somewhat smaller. Indeed, as we already 
mentioned, during the first oscillation the amplitude of 
the field $ is of order M p /20, which gives the Hubble 

constant H ~ ^/^f ir~ ~ 0.1m. Since dependence of 

the resonance condition on H is very weak (if 1 / 4 ), one 
may conclude that the regime of explosive reheating after 
inflation may occur only if the amplitude of oscillation 
satisfies the condition $ > m/g. Thus explosive decay 
ends at $ < m/g, i.e. at q < 1/4. 
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This means that preheating in this model cannot be- 
gin for $ < m/g, which would correspond to the narrow 
resonance regime. Narrow resonance may be important 
at the late stages of preheating jj], but at that stage 
one should take into account backreaction of the parti- 
cles produced at the previous stage of broad parametric 
resonance, so the theory of the narrow resonance at the 
end of preheating is much more complicated than the one 
contained in the previous subsection. 

In fact, efficient preheating often requires extremely 
large initial values of q. Indeed, the amplitude of the 
scalar field decreases during the expansion of the uni- 
verse much faster than H 1 / 4 , so for not very large initial 
values of q the condition ( |36| ) becomes violated before the 
resonance has enough time to transfer the energy of the 
oscillating field <fi into the energy of x-P&rticles. As we 
will show in Sec. IX. in the model under consideration 



preheating is efficient only if the initial value of q at the 
end of inflation is very large, qo > 10 3 . 

In the models with extremely large q the expansion of 
the universe makes preheating very peculiar: instead of a 
regular resonance process we encounter a rather unusual 
effect which we call stochastic resonance. 

Let us first look at the results of the numerical study 
of the development of broad resonance in an expand- 
ing universe, and try to interpret them. Note that at 
this stage we do not consider the effects of backreaction 
and rescattering of particles; we will discuss these effects 
later. Our main strategy here is to study a general pic- 
ture step by step, and then correct it later, because oth- 
erwise the physical interpretation of the processes which 
occur during preheating will remain obscure. 

First of all, let us consider Eq. for the mode Xk 
in an expanding universe with m 2 = 0, £ = in the 

asymptotic regime when a = (^-) 2 ^ 3 , and <&{i) — J^* mt ■ 
Strictly speaking, the last two conditions are satisfied 
only for sufficiently large t. However, if we begin counting 
time from the end of inflation, taking for definitcncss 
to = (which formally corresponds to the time after 
a quarter of one oscillation of the field <fi), then we will 
have an approximation which is sufficiently good for our 
purposes. With these definitions, the initial values of the 
field <j) and the parameter q in our calculations are given 
by 
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On the other hand, if one wants to investigate the situa- 
tion numerically, one can simply solve a combined system 
of equations for a(t), $(£) and Xk(t). We will not do it 
here because our main goal is to develop an analytical 
approach to the theory of preheating. 

The investigation of parametric resonance in an ex- 
panding universe can be simplified if instead of Xk one 
introduces the function Xu{t) = a 3 / 2 (t)xk(t), which is 
given by j~Xk{t) in our case. Then instead of ( pl| ) we 
have a much simpler equation 



where 
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(39) 



3 a. 
2 a 



£R. This term is usually 



very small. Indeed, we will consider here the case of light 
X particles, such that m x <C fc*, in which case one can 
simply neglect m x . Also, soon after the end of inflation 

one has H 2 = (f) 2 ^f < C ™ 2 - As a result, typically one 
can neglect the term A altogether. Eq. ( |38| ) describes an 
oscillator with a variable frequency u> 2 (t) due to the time- 
dependence of the background field <p(t) and a(t). As an 
initial condition one should take the positive-frequency 
solution X k {t) ~ e~ lUkt /y/2ujk~. 

The series of three figures in this section shows differ- 
ent stages of development of the fastest growing mode Xk 
in the broad resonance regime in an expanding universe 

2 ,2 

in the theory S for an initial value of the parameter 
g~3x 10 3 . Note that during the expansion of the uni- 
verse the amplitude of scalar field oscillations decreases 
approximately as t~ l . Therefore in order to illustrate the 
relative growth of the fluctuations of the field x with re- 
spect to the amplitude of the oscillating field 4> we show 
not the growing mode Xk itself, but its rescaled value 
Xk = Xk where to corresponds to the beginning of the 
calculation. One can construct an adiabatic invariant for 
Eq. ([38]), which has an interpretation of the comoving 
occupation number of particles n k in the mode k in an 
expanding universe: 



0Jk \X k \ 2 



\x k \ 2 



(40) 



Note that this function does not have any factors in- 
versely proportional to the volume a 3 . These factors will 
appear when we calculate the number density of particles 
in physical (not comoving) coordinates. 

In the beginning we have parametric resonance very 
similar to the one studied in the previous section, com- 
pare Fig. H and Fig. |]. As before, one can identify the 
periods when ^-particle production is most efficient with 
the intervals when the field ip becomes small. An impor- 
tant difference is that because of the gradual decrease in 
amplitude of the field (f> the effective mass of the field 
X and, correspondingly, the frequency of its oscillations 
decrease in time. As a result, in the beginning within 
each half of a period of oscillation of the field <fi the field 
Xk oscillates many times, but then it oscillates more and 
more slowly. 
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FIG. 4. Early stages of parametric resonance in the the- 
ory m £ in an expanding universe with scale factor a ~ i 2//3 
for g = 5 x 10 -4 , m = 10 _6 Af p . According to our conven- 
tions (0), initial value of the parameter q in this process is 
qo ~ 3 x 10 3 . Note that the number of particles in this pro- 
cess typically increases, but it may occasionally decrease as 
well. This is a distinctive feature of stochastic resonance in an 
expanding universe. A decrease in the number of particles is a 
purely quantum mechanical effect which would be impossible 
if these particles were in a state of thermal equilibrium. 

To understand the rather peculiar behavior of Xk and 
rik in this process, let us check in which resonance band 
our process develops. The number of the band in the 
theory of the Mathicu equation is given by n = \f~A~. In 
our case reheating occurs for A ~ 2q, i.e. n ~ %/2? ~ 
^? . Suppose we have an inflationary theory with m ~ 

10 _6 M p , and let us take as an example g ~ 10 _1 . Then 
after the first oscillation of the field, according to Eq. (jij), 
we have $(t) ~ M p /20, which corresponds to q ~ 
This gives the band number n~3x 10 3 . After another 
oscillation the amplitude of the field drops by a factor of 
two, and the band number decreases by a factor of two 
as well, down to n ~ 1.5 X 10 3 . 

In other words, even during a single oscillation the field 
does not remain in the same instability band of the Math- 
ieu equation. Instead of that it jumps over 10 3 different 
instability bands! The theory of a broad resonance in 
Minkowski space is much less explored than the theory 
of a narrow resonance, but the theory of a broad reso- 
nance in an expanding universe proves to be even more 
complicated. The standard method of investigation of 



resonance using the Mathieu equation in a single reso- 
nance band completely fails here. 

Still not everything is lost. Indeed, as we have found 
in the previous section, in the broad resonance regime 
particle production occurs only in a small vicinity of <f) = 
0, corresponding to integer and half- integer N — mtjl-K. 
Nothing depends on the exact way the field behaves 
at all other moments. In this sense the description of 
the process of particle production at 4> — is very robust 
with respect to change in the shape of the potential V{<j>) 
and of the equation describing the field \, insofar as it 
does not alter the behavior of either field at the stage 
when (pit) approaches zero. Therefore some (but not all) 
of the results related to the Mathieu equation can be 
useful for investigation of broad parametric resonance in 
an expanding universe even though the regime we are 
going to investigate is fundamentally different. 

One of the most important differences between broad 
resonance in Minkowski space and in an expanding uni- 
verse can be understood by inspecting the behavior of the 
phase of the functions Xk near the points where <f>(t) = 0. 
Indeed, Fig. || shows that near all points where = 
the phases of Xk are equal. The physical meaning of this 
effect is very simple: In order to open a swinging door 
by a small force one should apply it periodically, "in res- 
onance" with the motion of the door. 

However, in an expanding universe such a regime is 
impossible, not only because of the redshift of the mo- 
mentum -, but mainly because the frequency of oscilla- 
tions of the field Xk 1S proportional to $, which decreases 
in time. The frequency of oscillations of the modes Xk 
changes dramatically with each oscillation of the field 0. 
Therefore for large q the phases of the field Xk at suc- 
cessive moments when 0(i) = are practically uncorre- 
cted with each other. Using our analogy, one may say 
that the door is vibrating with a large and ever changing 
frequency, so it is very difficult to push it at a proper 
moment of time, and successfully repeat it many times 
in a row. That is why at some moments the amplitude 
of the field Xk decreases, see Fig. [|. 

This could suggest that broad parametric resonance 
in an expanding universe is simply impossible. Fortu- 
nately, this is not the case, for two main reasons. First 
of all, as we are going to show in the next section, even 
though the phases of the field Xfc at the moment when 
4>(t) = in an expanding universe with g> 1 are practi- 
cally unpredictable, in 75% of all events the amplitude of 
Xk grows at that time. Moreover, even if it were not the 
case, and the amplitude would grow only in 50% of all 
events, the total number of %-particles would still grow 
exponentially Indeed, as we will see, during each "cre- 
ative moment" <p(t) — in the broad resonance regime 
the number of particles at each mode may either decrease 
by a factor of O(10), or grow by a factor of O(10). Thus 
if we begin with 10 particles in each of the two modes, 
after the process we get 1 particle in the first mode and 
100 particles in the second. Therefore the total number 
of particles in this example grows by more than a factor 
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of 5. The theory of this effect is very similar to the theory 
of self-reproduction of an inflationary universe, where in 
most points the inflaton field rolls down, but those parts 
of the universe where it jumps up continue growing ex- 
ponentially Q. 

As a result, parametric resonance does take place. 
However, in order to describe it some new methods of in- 
vestigation of parametric resonance should be developed. 
We will do this in the next section. 

Stochastic resonance occurs only during the first part 
of the process, when the effective parameter q is very 
large and the resonance is very broad. Gradually the am- 
plitude of the field cf> decreases, which makes q smaller. 
Expansion of the universe slows down, the field stays in 
each resonance band for a longer time, and eventually the 
standard methods of investigation based on the Mathieu 
equat ion become useful again. As we will show in Sec. 
VII E , stochastic resonance ends and the standard meth- 

V4 



ods become useful after the first q /v27r oscillations, 
which may happen even before the effective parameter q 
decreases from qg ^> 1 to q ~ 1, see Eq. (|79|). One of 
the manifestations of the transition from the stochastic 
resonance to a regular one is a short plateau for Innfc 
which appears in Fig. | for 10 < t < 15. This plateau 
corresponds to the time when the resonance is no longer 
stochastic, and the mode Xk appears in the region of sta- 
bility, which divides the second and the first instability 
band of the Mathieu equation, see Fig. [?]. 




To get a better understanding of this effect one should 
continue our calculations for a longer period of time, see 
Fig. |[ At t > 15 the process does not look like a broad 
resonance anymore, but the amplitude still grows expo- 
nentially at a rather high rate until the amplitude of the 
field $ becomes smaller than m/g, which corresponds to 
q ~ 1/3 — 1/4. Soon after that the resonance ceases to 
exist and the amplitude stabilizes at some constant value. 



X, 



lnm 




FIG. 6. The same process during a longer time, which is 
shown in the units ^ , corresponding to the number of oscilla- 
tions N. The figures show the growth of the mode Xk. for the 
momentum k corresponding to the maximal speed of growth 
of rifc. In this particular case k ~ 4m. Towards the end of 
this period, after approximately 25 oscillations of the infla- 
ton field, the resonance ceases to exist, and the occupation 
number becomes constant. 




FIG. 5. The same process as in Fig. ^ during a longer 
period of time. The parameter q = 3 4 % decreases as t~ 2 
during this process, which gradually makes the broad reso- 
nance more and more narrow. As before, we show time t in 
units of — , which corresponds to the number of oscillations 
of the inflaton field. 



The time tf and the number of oscillations Nf at the 
end of parametric resonance in an expanding universe can 
be estimated by finding the moment when g$> as 4— | is 
equal to to: 
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gMp 
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(41) 



As one can check, this estimate for our case (to = 
IQ- 6 M P , g = 5 x 1CT 4 ) gives Nf ~ 26.5, which is in 
good agreement with the results of our computer calcu- 
lations shown in Fig. ||. A small disagreement (about 
10%) appears because our criterion for the end of the 
resonance <?$ ~ m was not quite precise: the resonance 
ends somewhat earlier, at g$ ~ 1.1m. 

This more exact result can be deduced from Fig. [7], 
which shows that the first instability band for k = 
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extends from q ~ 0.8 to q ~ 1/3. Therefore the growth 
of all modes with k <C m terminates not at g 2 ^ 2 /Am 2 ~ 
1/4, but slightly earlier, at g 2 $ 2 /4m 2 — 1/3. 

At the time t ~ t//2 one has q ~ 1. During the time 
from t f/2 to tf the resonance occurs in the first resonance 
band, the resonance is not very broad and there are no 
stochastic jumps from one resonance band to another. 
At the time just before tf/2 there was no resonance; the 
field was in the stability band between q = 1 and q = 2, 
see Fig. ^[ 
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FIG. 7. The structure of the resonance bands for the 
Mathieu equation along the line A — 2q, which correspond to 
excitations with k — in our model. The modes with small 
k are especially interesting because the momenta of the ex- 
citations are redshifted during the expansion of the universe. 
A small plateau at 10 < t < 15 on Fig. |B| corresponds to the 
time where stochastic resonance ceases to exist, all modes are 
redshifted to small k, and the system spends some time in the 
interval with 1 < q < 2, which is outside the instability zone. 
The last stage of the resonance shown in Fig. [] corresponds 
to the resonance in the first instability band with q < 1. 

An interesting effect which is shown in Fig. ^ is a slow 
growth of the amplitude which continues even after 
the resonance terminates and nfe becomes constant. This 
happens because the momentum of each mode gradually 
becomes smaller due to the expansion of the universe, 
and this leads to a growth of \k even though does not 
change. This is one of the examples which shows that 
in order to describe parametric resonance one should use 
proper variables such as rife, because otherwise one may 
get the incorrect idea that the resonance continues even 
for t > 25. 

If one ignores a small island of stability near t *~ 12, 
one may conclude that during the main part of the pro- 
cess the slope of the curve In nfe remains almost constant. 
In our case this corresponds to the exponential growth 
of the occupation number nfe with an effective parameter 
/ife ~ 0.13. This fact will be very useful for us later, when 
we will calculate the number of particles produced dur- 
ing the parametric resonance. Such a calculation is our 
main goal. It is also necessary in order to verify whether 
one should modify our resonance equations due to the 



presence of %-particles. As we will see, no modifications 
are needed for theories with g < 3 x 10~ 4 . However, 
for greater values of g (and in particular for the case of 
g ~ 5 x 10~ 4 discussed above) the resonance ends in a 



somewhat different way, see Sec. KB. 

In order to illustrate the stochastic nature of the res- 
onance in this theory, we will present here at sample of 
results for the resonance for several different values of 
the coupling constant g in the interval from 0.9 x 10~ 4 to 
10~ 3 . One might expect the results to change monoton- 
ically as g changes in this interval. However, this is not 
the case. The table contains the results concerning the 
initial momentum k (in units of m) corresponding to the 
fastest growing mode, the total increase of the number of 
particles In at the end of the resonance for this mode, 
the average value [i for this mode, and the time tf (the 
number of oscillations of the field 4>) at the end of the 
resonance: 
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k 


M 
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Inn/- 


0.9 x 10~ 4 


1.5 


0.1 


5 


6 


10~ 4 


2 


0.14 


5 


9 


1.1 x 10~ 4 


0.5 


0.17 


5.5 


12 


1.2 x 10~ 4 


1.5 


0.12 


6 


9 


1.3 x 10~ 4 


1 


0.13 


6.5 


11 


1.4 x 10~ 4 


2 


0.12 


7 


11 


1.5 x 10~ 4 


0.5 


0.18 


7 


17 


2 x 10~ 4 


3.5 


0.12 


11 


16 


3 x 10~ 4 


0.5 


0.14 


14 


27 


5 x 10~ 4 


4 


0.13 


24 


40 


10 _a 


6 


0.12 


48 


75 



Thus we see that the leading mode in this interval of 
the coupling constant has initial momentum comparable 
to m and slightly smaller than the typical initial width of 
the resonance fc*/2, which changes from 2m to about 5m 
for g changing from 10~ 4 to 10 -3 . The reason why k is 
usually (though not always) somewhat smaller than k*/2 
is very simple. The resonance is broad only during the 
first half of the time. Narrow parametric resonance which 
appears during the second part of preheating typically is 
more efficient for smaller k. We should note that for g > 
3 x 10 -4 , at the last stage of preheating one should take 
into account backreaction of produced particles, which 
make s the narrow resonance stage very short, see Sec. 
IX B. In such a case the resonance has the width fe*/2 in 
terms of the value of the momentum k at the beginning 
of preheating. 

Of course, investigation of the leading growing mode 
is insufficient: One should integrate over all modes with 
all possible k, which we are going to do later. However, 
the number of particles is exponentially sensitive to 
k. Therefore the main contribution to the integral will 
be given by the trajectories close to the leading one. It is 
similar to what happens, e.g., in the theory of tunneling, 
where one first finds the optimal trajectory corresponding 
to the minimum of action, and calculates e~ s along this 
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trajectory. Similarly, one can calculate the rate of growth 
of the total number of %-particles by finding the leading 
trajectory and calculating the average value of \x along 
the trajectory. 

The table clearly demonstrates that the effective values 
of /i and especially the final number of particles n k pro- 
duced by the resonance are extremely sensitive to even 
very small modifications of g, and change in a rather 
chaotic way even when g changes by only 10%. That is 
why we call this process "stochastic resonance." We see 
from the table that for g ~ 1CU 3 the occupation num- 
bers n k become incredibly large. It will be shown in Sec. 
[X that for g ~ 10~ 4 backreaction of created particles is 
not very important, but for g > 3 x 10~ 4 backreaction 
becomes crucial, because it does not allow the resonance 
to produce an indefinitely large number of particles. To 
investigate these issues we should first develop the the- 
ory of stochastic resonance, and then take into account 
backreaction. 



VII. ANALYTIC THEORY OF STOCHASTIC 
RESONANCE 

In this section we are going to develop a new method 
to study the time evolution of the eigenfunctions Xk(t) 
in the most interesting case of broad resonance. This 
method is based on the crucial observation made in the 
previous sections: In the broad resonance regime the evo- 
lution of the modes Xk (t) is adiabatic and the number of 
particles does not grow in the intervals when \<fi(t)\ > 0*. 
The number of particles changes only in the short inter- 
vals when \4>(t)\ < 0* < $. 

The quantum field theory of particle creation in a time 
varying background is naturally formulated in terms of 
adiabatic (semiclassical) eigenfunctions. This formalism 
is introduced in the next subsection. Then we will find 
the change of the particle number density from a sin- 
gle kick, when cf>(t) crosses zero at some time tj. For 
this purpose it is enough to consider the evolution of 
Xk(t) in the interval when <fi 2 (t) is very small, so it can 
be represented by its quadratic part oc (t — tj) 2 . This 
process looks like wave propagation in a time dependent 
parabolic potential. We can combine the action of the 
subsequent parabolic potentials to find the net effect of 
the particle creation. Using our formalism, we consider 
a toy model of broad resonance in Minkowski space, and 
broad resonance in an expanding universe, which turns 
out to have a stochastic nature. 



A. Adiabatic representation of the eigenfunctions 

The semiclassical, or adiabatic evolution of the eigen- 
function Xk(t) can be represented in a specific mathe- 
matical form. For this we adopt a physically transparent 



method to treat Eq. ( pq) for an arbitrary time depen- 
dence of the classical background field which was orig- 
inally developed by Zeldovich and Starobinsky |$7| for 
the problem of particle creation in a varying gravitational 
field. 

Let us represent solutions of Eq. ([38|) as products of its 
solution in the adiabatic approximation, exp (±i J dt u>), 
and some functions a(t) and (3(t): 



aV 2 Xk {t)=X k (t) 
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(42) 



An additional condition on the functions a and (3 can be 
imposed by taking the derivative of Eq. ( ^2| ) as if a and 
/3 were time- independent. Then Eq. ( j42] ) is a solution of 
Eq. (|38|) if the functions a k ,/3 k satisfy the equations 
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In terms of classical waves of the x-held, quantum effects 
occur due to departure from the initial positive-frequency 
solution, therefore the initial conditions at t — > are 



1, (3 k = 0. Normalization gives \a k \ 2 — \P k \ 2 = 1. 



The coefficients a k (t) and f3 k (t) in our case coincide 
with the coefficients of the Bogoliubov transformation of 
the creation and annihilation operators, which diagonal- 
izes the Hamiltonian of the x -ne ld at each moment of 
time t. The particle occupation number is n k = \Pk\ 2 , 
see Eq. ( [io"|) . The vacuum expectation value for the par- 
ticle number density per comoving volume is 
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(44) 



n x ) in the 



In this section we will calculate (3 k , n k and 
non-perturbative regime of broad resonance, where all of 
these values can be very large. 

It is instructive to return in the framework of this for- 
malism to the simpler perturbative regime which we dis- 
cussed earlier in Sec. III. Assuming \(3 k \ <C 1, from Eqs 
([43]) one can obtain an iterative solution: 



t' 



dt — exp (—2i 
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(45) 



Using uj(t) = y (^) 2 + g 2 $ 2 sin 2 mt, we can evaluate Eq. 

( [45| ) containing an oscillating integrand by the method 
of stationary phase . In the case of the massive scalar 
field decaying via the interaction g 2 a(f>x 2 , the dominant 
contribution is given by the integration near the moment 
ife, where a(t k ) = As we already mentioned, this cor- 
responds to the creation of a pair of massless x-particles 
with momentum k = ■^a(t k )m from an inflaton with mass 
(energy) m at the instant t k of the resonance between 
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the mode k and the background field. The decay rate 
of the inflaton field calculated with this method can be 
described by Eq. ( |l2| ) . 

For the interaction (jxp — > xXi the process in the regime 
|/3fc| <C 1 can be interpreted as creation of a pair of \~ 
particles with momentum k = a(tk)m from a pair of 
massive inflatons with energy m each. The decay rate 
of the massive inflaton field in this case rapidly decreases 
with the expansion of the universe as ^^(a 4 Px) ^ a ~ 6 - 
Therefore a complete decay of the massive inflaton field in 
the theory with the </> 2 x 2 -interaction is impossible. One 
should have additional terms such as g 2 a4>x 2 or hifnptfi. 
This is a very important conclusion which we already 
discussed in Sec. 



[II. 



B. Interpretation of parametric resonance in terms 
of successive scattering on parabolic potentials 

We suggest a new analytic method to solve approxi- 
mately the basic equations (|2l|) and ( |38| ) for the eigen- 
functions Xk which correspond to the x-particles created 
by the oscillating inflaton field <p(t). This method is 
rather general; it can be applied to many models of pre- 
heating. One may also apply it to the idealized case 
when the universe does not expand and backreaction is 
not taken into account. In the cases where the equation 
for the modes Xk can be reduced to an equation with pe- 
riodic coefficients (including the Mathieu equation), our 
method accurately reproduces the solution of this equa- 
tion, and gives us an interesting insight into the physics 
of parametric resonance. This method is rather powerful; 
it enables one to investigate some features of the regime 
of broad parametric resonance which, to the best of our 
understanding, have not been known before. 

In the realistic situation which we study in this paper, 
when the expansion of the universe as well as the backre- 
action are taken into account, in some models (e.g. non- 
conformal theory) the equation for the modes Xk cannot 
be considered as an equation with periodic coefficients, 
and the analysis based on standard stability/instability 
charts is not applicable. This is the situation where our 
method will be especially useful. 

Let us consider the general equation (|3|). As we no- 
ticed, the eigenfunction Xk(t) has adiabatic evolution be- 
tween the moments tj, j — 1,2,3, where the inflaton 
field is equal to zero 4>{tj) = 0, (i.e. twice within a pe- 
riod of inflaton oscillation). The non-adiabatic changes 
of Xk (t) occur only in the vicinity of tj. Therefore we 
expect that the semiclassical solution ( |42"|) of Eq. (|3q ) 
is valid everywhere but around tj. Let the wave Xk(t) 
have the form of the adiabatic solution ([l2|) before the 
scattering at the point tj 



the coefficients oci and (3j, are constant for tj—\ 



(46) 



Then after the scattering, Xk(t), within the interval tj < 
t < ij+i, has the form 



jiii 
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and the coefficients and f3i +1 are constant for tj < 



t < t 



j'+i 



Eqs. (46) and ( |47|) are essentially the asymptotic ex- 
pressions for the incoming waves (for t < tj) and for the 
outgoing waves (for t > tj ), scattered at the moment 
tj. Therefore the outgoing amplitudes cP k +1 , /3 J k +1 can be 
expressed through the incoming amplitudes oA, 01 with 
help of the reflection and transmission Dk amplitudes 
of scattering at tj : 



j'+i - 

k e 




(48) 



Here 9° k = J dt oj(t) is the phase accumulated by the 
o 

moment tj. 

Now we specify the scattering at the moment tj. The 
interaction term g 2 (f> 2 (t) in Eq. ( p8| ) has a parabolic form 
around all the points tj'. g 2 4> 2 (t)~ g 2 § 2 m 2 (t — tj) 2 = 
k*(t — tj) 2 , where the current amplitude of the fluctu- 
ations $ is defined in (0), and the characteristic mo- 
mentum = yjg^m. In the general case depends 
on time via the time dependence of $ cx a~ 3 / 2 . Figure 
8 illustrates two possible outcomes of the scattering of 
the wave Xk(t) on the parabolic potential near zeros of 
the function g 2 (f> 2 (t). Depending on the phase of the in- 
coming wave, the corresponding number of particles may 
either decrease or grow. 
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FIG. 8. The change of the comoving particle number nu 
due to scattering at the parabolic potential, calculated from 
Eq. (|5o|). The dotted lines show the sequence of the parabolic 
potentials g 2 cj) 2 (t) w g 2 & 2 m 2 {t — tj) 2 where scattering occurs. 
Time is given in units of — . The number of particles can 
either increase or decrease at the scattering, depending on 
the phase of the incoming wave 
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First, let us consider the mode equation around a sin- 
gle parabolic potential. In the vicinity of tj the general 
equation (p8[) is transformed to the equation 



d 2 X k (k 2 
dt 2 



+ [■-+gH 2 m 2 (t-t ] f 



X k = 



(49) 



For simplicity we introduce a new time variable r = 
fe* (t — tj ) and a scaled momentum k = -|- . Notice that 

k' 2 = (A k — 2q)/2-Jq. In general, fc* and k depend on 
tj through a(tj), and should be marked by the index j, 
which we drop for the moment. Then Eq. ( f49| ) for each 
j is reduced to the simple equation 

d 2 X 

+ {n 2 + t 2 ) X k = . (50) 

The asymptote of this equation, which corresponds to 
the incoming wave, is matches to the form ([46|). The 
asymptote corresponding to the outgoing wave matches 
the form (p7|) . Therefore the reflection R k and trans- 
mission D k amplitudes of scattering at tj are essentially 
the reflection and transmission amplitudes of scattering 
at the parabolic potential. Thus the problem is reduced 
to the well-known problem of wave scattering at a (neg- 
ative) parabolic potential |3S| ], which we consider in the 
next subsection. 



C. Particle creation by parabolic potentials 

A general analytic solution of Eq. (^0|) is the lin- 
ear combination of the parabolic cylinder functions |}9| : 

W ±y/2rj . The reflection R k and transmission 

D k amplitudes for scattering on the parabolic potential 
can be found from these analytic solutions: 



Rk — — 



D k 



VT 



-up k 



VT 



(51) 



(52) 



(53) 



The angle if depends on the momentum k. Notice the fol- 
lowing properties of these coefficients: R k = —iD k e~ 



where the angle tp k is 
(fik = argT 



1 + in 2 \ k 2 f , 2 
1 + hi — 



l-Rfel + \ D k\ = !• Substituting Q51p and (|52|) into (|48|), 
we can obtain the evolution of oci, 01 amplitudes from a 
single parabolic scattering in terms of the parameters of 
the parabolic potential and the phase 9 k only. 
The mapping of crjj,, (3 J k into crjj, +1 , P 3 k +1 reads as 



(Y 



3+1 



(54) 



VI + e- 



le 



— S k +2i0 



VT 



Since the number density of x-particles with momen- 
tum k is equal to n k — \(3k(t)\ 2 , from Eq. ( |54] ) one 
can calculate the number density of outgoing particles 
= l/3fc +1 | 2 a ft cr the scattering on the parabolic po- 
tential out of n 3 k = \(3 k \ 2 incoming particles: 



,3+1 



l + 2e" 



2e"2 



« Vl + e-« 2 n{ (1 + n{) sin 6{ ot , (55) 



where the phase 8 3 tot = 29 k — ip k + arg/3^, — arg af.. 

Before we apply the formalism (Q) and (g5|) to spe- 
cific models, we shall analyze these generic equations. 
Although we did not specify yet the phase 9l ot , we al- 
ready can learn a lot from the form (|55| ) . First of all, the 
number of created particles is a step- like function of time. 
The value of n k is a constant between two successive scat- 
terings at points tj and tj+i. The number of particles is 
changed exactly at the instances tj in a step-like man- 
ner, in full agreement with the exact numerical solution, 
see Figure 4. The effect of particle creation is signifi- 
cant if ttk < 1, otherwise the exponential term e 7rK 
suppresses the effect of particle accumulation. This gives 
us the important general criterion for the width of the 
resonance band Q : 



.4 



2V? 



(56) 



where A = 



2<7, q 



4m 2 



Equivalently, one can 



write this condition in the following form: 



k 2 

- 2 < k 2 /n 



(57) 



This estimate of the resonance width k < fc* /Vtt practi- 
cally coincides with the estimate k < k*/2 ( p4| ) derived 
in Sec. [y] by elementary methods. 

Next, let us consider the large occupation number 
limit, n k 3> 1. From Eq. (|5q) we derive 



,3 + 1 



l + 2e" 



2sinCt 



VT 



(58) 



The factor in the r.h.s. of this equation depends on 
the coupling constant g through k 2 oc g~ . This de- 
pendence has the structure exp (— 1/g), which is a non- 
analytic function of g at g = 0. Therefore the number of 
particles generated in the broad resonance regime cannot 
be derived using a perturbative series with respect to cou- 
pling parameter g. Thus formula ( pq ) clearly manifests 
the non-perturbative nature of the resonance effects. 
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The growth index /j,k is defined by the formula 

n l +1 ="fcexp(27r/4) . (59) 
Comparing ((5^) and ([58]) we find 

H j k = — In + 2e _7rK2 - 2 sin^ ot e - '* 2 \/l + e-™ 2 ) . 

(60) 

The first two terms in Eq. (|58f) correspond to the effect 
of spontaneous particle creation, which always increases 
the number of particles. The last term corresponds to 
induced particle creation, which can either increase or 
decrease the number of particles. At first glance it looks 
paradoxical that the number of particles n k created from 
the time-varying external field can not only increase but 
sometimes decrease, i.e. the growth index \i k can be not 
only positive but sometimes negative. Indeed, it is well 
known that if the |m)-state of the quantum field x corre- 
sponds to n particles, then the number of particles in the 
|oui)-state due to the interaction with the external field 
will always be greater than n. This is how to resolve the 
paradox: the particles created from the vacuum by the 
time-varying external field are not in the n-particle \in}- 
state but are in the squeezed |m)-state. In this case the 
interference of the wave functions can lead to a decrease 
of the particle number. 

The whole effect of the particle production crucially 
depends on the interference of the wave functions, i.e. 
the phase correlation/anticorrelation between successive 
scatterings at the parabolic potentials. The maximal 
value of ji is reached for positive interference when 
sin# tot = —1 and is equal to = — In (l + %/2) ~ 0.28, 
see also [|o|. The typical value of fj, corresponds 
to sin#i 0t = and is equal to = ^hi3 w 0.175. 
The value of /j, is negative for negative interference when 
sin6> tot = 1. Therefore the behavior of the resonance 
essentially depends on the behavior of the phase 9 k as 
a function of k for different time intervals j, see Fig. 
8. In the case of a fixed amplitude of the background 
field <I>(t) = const and a(t) — const, the phases 9 k do 
not depend on time but only on k. In this case we ex- 
pect the existence of separate stability and instability 
/c-bands. However, this separation is washed out as soon 
as the phases 9 k are significantly varying with time due 
to changes in the parameters of the background field, for 
instance, in $(i) and a(t). 

Now we estimate the net effect of particle creation after 
a number of oscillation of the inflaton field. Eqs. (p8|) and 
([55|) are recurrence relations for the a J k and /3 3 k coefficients 
and for the number of particles n k after successive actions 
of the parabolic potentials centered at ti , t%, .... To find 
the number of particles created up to the moment tj , one 
has to repeat the formulas j-times for the initial values 
a k = 1, P J k = 0, n k = and a random initial phase 9 k . 

After a number of inflaton oscillations, the occupation 
number of x-particles is 



n k (t) 



(61) 



where we convert the sum over (ir k to an integral over 
/ifc(i). In some cases the index u k (t) does not de- 
pend on time. In a more general case one can replace 
Mfe(^) by an effective index n° k s defined by the relation 
J dt/Jk(t) = which, for brevity, we will write sim- 

ply as u k t. Then the total number density of created 
particles is given by 



n x (t) 



(2na) 3 



d 3 k n k (t) = 



47r 2 a 3 



dk k 2 e 2m ^ kt . 

(62) 



The function fi k has a maximum ^ max = /j, at some k 
■ The integral ( |62| ) can be evaluated by the steepest 
descent method: 



n x (t) 



1 



1 



47r 2 a 3 



2tt mt ji k 



8n 2 a 3 y/Ttfimt 



(63) 



where n k is the second derivative of the function [i k at 
k = k m which we estimated as fi k ~ 2/^/Afc 2 , Ak be- 
ing the width of the resonance band. Thus the effect of 
particle creation is defined by the leading value of the 
growth index /x, by the leading momentum k m and by 
the width of the resonance band Afc. In practice typi- 
cally k m ~ Ak ~ fc*/2, so we can use an estimate 



n x (t) 



647r 2 a 3 y / 7r7u77i 



(64) 



In order to calculate n(t) one should find the values of 
the parameters /i and fe*. 

In what follows in this section we will apply the general 
formalism of successive parabolic potentials first to the 
toy model without the expansion of the universe, where 
a(t) = const and $(t) = const, in the case of broad res- 
onance, 5 > 1. We will find the resonance zones and 
the number of particles which would be created in such a 
model. Then we consider a realistic case with the expan- 
sion of the universe taken into account. It turns out that 
the resonance in an expanding universe is very different 
from that without expansion. 



D. Broad parametric resonance without expansion 
of the universe 



Let us apply the general formalism of the previous sub- 
section to the toy model neglecting the expansion of the 
universe. This is equivalent to taking a(t) = 1. Thus, 
we will study the evolution of the eigenfunctions in the 
case with fixed values of the background parameters and 
without backreaction of created particles. In this case 
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Eq. ( |2"l| ) is reduced to the standard Mathieu equation 



( |23| ) with Ak 



2q, i 



— 9_ 



z = mt. 



As we saw in Sec. VI, for the realistic situation with 



the expansion of the universe the Mathieu equation is 
applicable only at the last stages of the resonance when 
q < 1. For 9 > 1 this equation has only a heuristic 
meaning for our problem. 

For the Mathieu equation with a large value of q (which 
is a constant in this subsection) we have the broad reso- 
nance regime. In this case the parameters k 2 and ipk of 
matrix (|54| ) are time-independent, i.e. they are the same 
for different j. The phase di is simple: di = Ok' 3- Here 
8k = \J dt uik is the phase accumulating between two 
successive zeros of <j>(t), i.e. within one half of a period of 
the inflaton oscillations, ir/rn, so that 9k = J m dt u>k- To 
find erf, and j3 k we have to apply the same matrix ( p34| ) 
j times. We are mainly interested in the regime with 
a large number of created particles, n k = \f3 k \ 2 3> 1. In 
this regime \a k \ ~ \(3 k \,so a k and (3 k are distinguished by 
their phases only. In this case there is a simple solution 
of the matrix Eq. (p4|) for an arbitrary j: 



1 



V2 



It) 



,(irfi k +i8 k )-j 



,{-Kiik-i6h)-i 



(65) 



(66) 



where 1? is a constant phase. In principle, it is possible to 
construct not only the asymptotic solution (|65|), (|66|), but 
the general solution which starts with (3 k = 0. However, 
the general solution very quickly converges to the sim- 
ple solution (^35|), (^), which contains all the physically 
relevant information. From (|66| ) the number of particles 
created by the time t « — is 



_ __ e 2irfj. k -j _ }_ e 2fj, k mt 

2 2 



(67) 



where \ik from (pq), (pq ) is indeed the growth index. 
Substituting the solution ( |65| ) , (|6^) into Eq. (|54j) , we get 
a complex equation for the parameters [ik and 9k 



(68) 



Alongside the solution ( |65| ) , ( pq ) , there is another asymp- 
totic solution of the matrix equation (|54|) 



1 



V2 



,(irii k +iQ k +iir)-j 



(69) 



Excluding the phase i9 from the complex equations ( pq ) 
and (0), it is easy to find a single equation for the growth 
index /z/. valid for both solutions :[] 



s*"* =|cos(fl fc - V j fc )|vl + e = 



(72) 



+ ^/(l + e"™ 2 ) cos 2 (0 fe - p fc ) - 1 • 



In the instability bands, the parameter \ik in Eq. (|72 
should be real. Therefore, the condition for the momen- 
tum fc to be in the resonance band is cos(9k — fk) > 
or 



tan(0 fc - cpk)\ < e~ 



(73) 



To further analyze the constraints on the width ( |73] ) 
and strength ([72]) of the resonance, we should find how 
the phases 9k and (fk depend on the momentum fc. The 
angle <fk as a function of fc is defined by Eq. (|53|). For 
the phase 9k we have 



dt^/W+g^if) 



2o$ k 2 / o$ 
In - 



2 \ IflK' 

fc2 



41n2 + 1 



(74) 



A^/qr, 



\q 



41n2 + 1 



To obtain these estimates we used the condition that 
ft2 ^ f° r * ne resonan t modes. In Eq. ( f74| ) we pre- 
sented 9k in two equivalent forms: first in terms of the 
physical parameters g, $, and k, and second in terms of 
the parameters q and fc. Combining Eqs. (f74|) and j5^ ) 
for the phases and fk, we can find how 9k — fk de- 
pends on fc. The leading term in 9k — fk for large values 
of q is the term — ^\fq which does not depend on fc. 
Substituting 9k — fk into Eq. ( [73] ) we get the equation 
for the width of the resonance explicitly in terms of k 
for a given parameter q. Eq. ([73]) transparently shows 
the presence of a sequence of stability/instability bands 
as a function of k. Typical half- width of a resonance 
band is fc 2 ~ O.Ifc 2 . Substituting 9 k — fk into Eq. ^7%), 
we find the strength of the resonance as a function of fc. 
The effect of amplification is not a monotonic function 
of q. The strongest amplification is realized for discreet 
values of the parameter q: q — (^j-) , where n is an 
integer. For this case [ik has a maximum at k = 0. We 
can illustrate our results graphically for this case, since 



aj _ J_.il) . Jirp, k -i9 k -iTr)-j 

Pk - 7f e e 



with the condition 



(-Kli k +i9 k 



t Notice that the number of particles calculated with Eq. 

(70) (Q) is in an agreement with the general formula jHq 1 ). From 
the definition of k ,tot and the solutions a 3 k and /3 k we 
have k tot ~ fk — Therefore from the complex equa- 
ti ons (pq 1 ) and ( |7l] ) we have additionally that cos(8 k ,tot) = 

(71) + sin(6»fc-v5 fc )- 
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the function argr( 1+ " ) involved in the expression for 
ifk (53) has a particularly simple form for k 2 « 1: 



1 



argT 
Then we have 



-0.982 k 



■(\nq + 9.474) 



(75) 



(76) 



The function [ik derived with the formulas (jT^) and ( |76| ) 
is plotted in Fig. 9 for q — (327r) 2 . We also plot u,k de- 
rived numerically from the Mathieu equation (|30|). We 
conclude that the predictions of the analytic theory de- 
veloped here for the Mathieu equation with large q are 
rather accurate. 
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FIG. 9. The characteristic exponent [i^ of the Mathieu 

The 



equation (^) as a function of k 2 



4, for q = (32tt) 2 
dotted curve is obtained from a numerical solution. Two in- 
stability bands are shown. The solid curve for these instability 
bands was derived analytically with Eqs. ( fri| ) and ( |76| ) where 
the simple approximation ( |75| ) was used. The numerical and 
analytical results are in a perfect agreement for the first band 
where the approximation (75) is accurate. By improving ex- 
pansion (75), one can reach similar agreement for the higher 
bands as well. 



E. Stochastic resonance in an expanding universe 

Let us consider the creation of x-particles by harmonic 
oscillations of the inflaton field in an expanding universe. 
Due to the expansion of the universe, there are few com- 
plications in Eq. ( pl| ) for the modes Xfc m an expand- 
ing universe in comparison with the Mathieu equation. 
The effect of the term 3Hy can be eliminated by using 
Xk = a 3 / 2 Xfc, see Eq. (p8[). The redshift of momenta 
k — > should be taken into account, especially at the 
latest stages. The most important change is the time- 



dependence of the parameter q 



4m 2 



q cx t 2 (x N 



For the broad resonance case where q ^> 1, this parame- 
ter significantly varies within a few inflaton oscillations; 
hence, the concept of the static stability/instability chart 
of the Mathieu equation cannot be utilized in this impor- 
tant case. 

Surprisingly, the most interesting case when the pa- 
rameter q is large and time-varying can also be treated 
analytically by the method of successive parabolic scat- 
terings. Indeed, the matrix mapping for the a k and f3 k 
developed in subsections B and C is also valid in the case 
of an expanding universe. Let us consider the phase ac- 
cumulating between two successive zeros of the inflaton 
field: 



dt 



rn 

gMp 

5 m j 
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T 
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In 



g<L> 



41n2 + 1 



(77) 



0{k 2 ) 



where we used Eq. (|) for the amplitude of oscillations, 
<&, as a function of the number of oscillations, j « 2N. 
If the initial value \q^i = yfqo is large, then variation of 
the phase 58 k between successive scatterings due to the 

, or in terms of the number 



j-dependence is S9 k 
of oscillations 



gM p 



"" K 20mN 2 2N 2 ' ^ 
The crucial observation is the following: for large initial 
values of q, the phase variation 89k is much larger than 
7r for all relevant k. Therefore, all the phases Qi in Eqs. 
(Jl8|) and ( [So] ) in this case can be considered to be random 
numbers. For given qo, the phases are random for the first 



-/Vstoch 



1/4 



2tt 



(79) 



oscillations. For example, for qQ = 10 6 the phases are 
random for the first dozen oscillations, and for qo = 10 9 , 
neglecting backreaction effects, the phases would be ran- 
dom for the first hundred oscillations. During this time 
each mode experiences chaotic behavior in the standard 
terms of the theory of chaotic systems [Q: a small 
change in the values of parameters and/or initial con- 
ditions can lead to large chang es in the final results. 



We will show in Sec. IX B that the backreaction of 



created particles leads to an exponentially rapid decrease 
of q down to q ~ 1/4 at the last moments of preheating. 
This means that the parameter q in this regime remains 
very large and phases remain random until the very last 
stages of preheating. 

The stochastic character of the phases, 9 k , significantly 
simplifies the analysis of the matrix equation (|48|). In- 
deed, since there is no memory of the phases, each map- 
ping can be considered as independent of the previous 
ones. 
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As we see in Eq. ([38|), the number of created particles 
depends on the phase 9 to t = (pk + 29l + arg (3k — arg a% . In 
principle, from the matrix equation (E ) one can derive a 
series of equations which allow one to express the phases 
a,rgf3k, argctfc, and eventually 9 tot through the random 
phase 9 k . 

For qualitative analysis we simply assume that 9 to t is a 
random phase. As a result the number of particles n^ k +1 
obeys the recursion equation 



1 + 2e" 



2sin( 



1 



(80) 



where 9 is a random phase in the interval (0, 2ir), and k 2 
is slowly changing with j as k 2 = oc j _1/3 . 

Eq. ( p0| ) defines the number of particles at an arbitrary 
moment as a function of the random phase. Therefore, n k 
is a random variable which can either increase or decrease 
depending on the realization of the phase. Qualitatively, 
each mapping corresponds to one of the two possibilities 
depicted in Fig. 8. Therefore, the whole process of parti- 
cle creation is the superposition of elementary processes 
where nk jumps up or down. This explains the random 
behavior of nk in Fig. 5. On average the number of par- 
ticles is amplified with time, i.e. nk increases more often 
than it decreases. 

Stochastic resonance is different in many aspects from 
the usual broad parametric resonance of the Mathieu 
equation, considered in the previous subsection. Let us 
investigate the basic features of the stochastic resonance. 
First, the structure of Eq. (|8(]) does not imply the exis- 
tence of separate stability or instability bands. Indeed, 
the loss of the phase interference appears for any k within 
the broad interval k < fc* , where the coefficients of the 
mapping (|8^) are not exponentially suppressed. There- 
fore, as one can see by comparison of Figs. || and [h], the 
stochastic resonance is significantly broader (almost by 
an order of magnitude) than each of the stability zones 
of the Mathieu equation, Afc ~ fe*. It makes stochas- 
tic resonance more stable with respect to possible mech- 
anisms which, in principle, could terminate parametric 
resonance. For instance, the conclusion that the g 2 (j> 2 x 2 
interaction can terminate broad parametric resonance in 
Minkowski space-time |29| cannot be easily generalized 
to the case of an expanding universe, where the broad 
resonance is stochastic and much wider. 

Second, the exponent [ik is also a random variable: 



/4 = ^lnM 



2e~ 



- 2sin( 



1 



(81) 



The functional form of \ik for stochastic resonance is 
different from that for broad parametric resonance. It 
changes with every half period of the inflaton oscilla- 
tions. An example of \ik calculated at intermediate stage 



of stochastic resonance (for j = 10) with the initial value 
of the parameter q ss 10 4 is plotted in Fig. [l(]. 




FIG. 10. The characteristic exponent \i k of the mode Eq. 
( p8[ ) in an expanding universe as a function of n 2 for the ini- 
tial value of the parameter q = (327r) 2 ~ 10 4 , obtained from 
a numerical solution. The curve is obtained at the time after 
the first 5 oscillations, which corresponds to fi k with j — 10. 
The envelope of the curve is obtained from Eq. ( |8l| ) by taking 
there sin# = ±1. We see that there is a complete agreement 
between the analytical prediction of the amplitude of ^ k ( Jsi| ) 
and the results of the numerical investigation. Contrary to 
the static case of Fig. [], the resonance is much broader, there 
are no distinguished stability/instability bands, and for cer- 
tain values of momenta the function n k is negative. During 
the stochastic resonance regime, this function changes dra- 
matically with every half period of the inflaton oscillations. 
Comparison of Figs. ^ and |Io| shows that it is incorrect to 
use the structure of the resonance bands of the static Mathieu 
equation for investigation of the stage of stochastic resonance, 
unless one is only looking for a very rough estimate of /i. 



Equation ( |3l| ) implies that for ttk 2 <g; 1 the value of 
fi k is positive (i.e. the number of particles grows) for 
f < 9 < . This occurs for one quarter of all possi- 



ble values of 9, in the range of — tt < 9 < n. Therefore, 
positive and negative occurrences of \ik for k <C 1/ \pK ap- 
pear in the proportion 3:1, so that the probability for the 
number of particles to increase is three times higher than 
the probability of its decreasing; see Sec. VI. Computer 
simulations of this process confirm this result. However, 
there will be also a "natural selection effect" : among all 
modes \k there will be some modes for which positive oc- 
currences of [ik appear more often than in the proportion 
3:1, and these modes will give the dominant contribution 
to the total number of produced particles. The typical 
mean value of the characteristic exponent is [ik ~ 0.13, 
but the actual number is very sensitive to even a very 
small change of parameters; see the table in Sec. VI. 



Based on the central limit theorem, we expect that the 
statistics of the random variable nk obey the log-normal 
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distribution in the regime of the stochastic resonance. 

From Eq. (^lj) one could expect that the suppression 
of particle production occurs not at k 2 > tt^ 1 , but at 
k 2 > 2ir~ 1 . However, the situation is more complicated. 
As soon as the second term under the logarithm becomes 
small, the probability for the number of particles to in- 
crease becomes equal to the probability of its decreasing, 
so the process of particle production becomes much less 
efficient. 

The stochastic resonance occurs for N st0 ch oscillations 
of the inflaton field defined by Eq. (f79f). When the pa- 
rameter q decreases because of the expansion of the uni- 
verse and becomes smaller than 0(1), which happens for 
N > iVgtoch, the resonance becomes very similar to the 
usual parametric resonance with q < 1. However, at some 
stage it may become necessary to correct this description 
by taking into account backreaction of the created parti- 
cles. 



VIII. RESONANCE, BACKREACTION AND 
RESCATTERING 

Until now we have treated the field x as a test field 
in the presence of the background fields 4>{t) and a(t) 
which have independent dynamics. We found the effect 
of the resonant amplification of Xfe(i), which corresponds 
to the exponentially fast creation of n x particles. As we 
have seen, the resonance in an expanding universe in the 
beginning may be very broad, then it becomes narrow, 
and then eventually disappears. 

Because of the exponential instability of the x field, we 
expect its backreaction on the background dynamics to 
gradually accumulate until it affects the process of reso- 
nance itself. Therefore the development of resonance is 
divided into two stages. At the first stage of the process, 
the backreaction of the created particles can be neglected. 
As we will see, this stage is in fact rather long, and if the 
initial value of q was small enough (q < 10 3 ) preheating 
may end before the backreaction becomes important (see 
also [p8[). However, if qo is greater than about 10 3 , then 
at some moment the description of the parametric reso- 
nance changes. We enter the second stage of preheating 
where the backreaction should be taken into account. In 
what follows we will treat the first and second stages of 
preheating separately. 

There are several ways in which backreaction can alter 
the process. First of all, interaction with particles created 
by parametric resonance may change the effective masses 
of all particles and the frequency of oscillation of the 
inflaton field. Also, scattering of the particles off each 
other and their interaction with the oscillating field (j>(t) 
(we will vaguely call both processes "rescattering" ) may 
lead to additional particle production and to the removal 
of previously produced particles from the resonance. 

In our model there will be two especially important 
effects. First, %-particles may change the frequency m 



of oscillations of the field 4>{t). This may increase the 
value of m in the mode equation, which can make the 
resonance narrow and eventually shut it down. 

The second effect is the production of ^-particles, 
which occurs due to interaction of x-particles with the os- 
cillating field 4>{t). One can visualize this process as scat- 
tering of x-particles on the oscillating field cf>(t) . In each 
act of interaction, each x-particle takes one ^-particle 
away from the homogeneous oscillating field 4>(t). When 
many 0-particles are produced, they may change the ef- 
fective mass of the field x, making ^-particles so heavy 
that they no longer can be produced. Also, scattering, 
when it occurs for a sufficiently long time, can destroy 
the oscillating field 4>{t) by decomposing it into separate 
^-particles. 

In this section we will derive the general set of equa- 
tions which describe the self-consistent dynamics of the 
classical homogeneous inflaton field cf>(t) , as well as the 
fluctuations of the fields x an d 4>- We will then discuss 
different feedbacks of the amplified fluctuations. In par- 
ticular, we will check the energy balance between the 
background homogeneous inflaton field <fi(t), the fluctua- 
tions x(i; x ); and the fluctuations 0(i, x). 

A. Self-consistent evolution of (f> and \ fields 

We can describe all of these effects within a full set 
of self-consistent equations. The Fricdmann equation for 
a universe containing classical field <f>(t) and particles x 
and <j) with densities p x and is 

■AH 2 = ^(\4> 2 + \m 2 4> 2 + Px + P(k ) , (82) 

where p x and p^ are the energy densities of x-particles 
and </>-particles respectively. 

The mode Eq. (||) for X k {t) = a 3 / 2 (t)xk(t) now 
should include a term describing the coupling between 
X and <fi fluctuations: 

X k (t) + + ,g 2 $ 2 sin 2 mt\ X k {t) = 

- J dt'X k (t')IL x (t,t';k) , (83) 

where the polarization operator for the field Xk — 
a- 3 ' 2 X k is Il x (t,t';k) = J d 3 j;e ,k ( x - x '>n x (t,f';x-x'). 

We will also consider quantum fluctuations of the in- 
flaton field S(f>(t, x) = 0(t,x) — <fi(t) which can exist on 
top of the homogeneous inflaton condensate <fi(t). The 
mode equation for tpk(t) = a 3 ^ 2 Sc/)k(t) is 

$k{t) +(^2+ ™ 2 ) ¥>*(*) = - J dt'ipkit')^, t'; k) , 

(84) 
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where H<j,(t, t';k) is a corresponding polarization opera- 
tor for the field 8<fik(t) = a~ 3 / 2 (fk(t). The equation for 
the homogeneous condensate <j)(t) is 

4>(t) + 3Hj>(t) + m 2 ^{t) = -r^(t) = -nj(t) . (85) 

Here r^(i) is the tadpole diagram, representing the 
derivative of the effective action of the field <f> (not the 
decay rate!). The one-loop diagram representing T^(t) 
is shown in Fig. The thick line corresponds to the 
exact solution of the classical equation of motion of the 
field x m the external field <fr. 




4>(t) 

a) b) 

FIG. 11. The one-loop diagrams for the contribution of 
X-particles to the effective action of the field (j> (Fig |ll|a) and 
to its derivative T^,(t) (Fig |ll|b). The thick line corresponds 
to the Green function of the ^-particles in the external field 

To get an expression for the polarization operator of 
the field cf>, one should differentiate the effective action 
twice with respect to the scalar field Scf>. The result can be 
represented as a sum of two polarization operators shown 
in Fig. |l2|. U\ can be identified with the contribution of 
the fluctuations of the field x to the mass squared of the 
field <fi: Am 2 = g 2 (x 2 )- Note that it is directly related 
to r^: r0 — — II 1 4>. The polarization operator 11^ has a 
more complicated structure; it contains an external scalar 
field 4>(t) in each of its vertices due to the interaction 
g 2 8q>cj ) (t)x 2 . 

The self-consistent dynamics described by Eqs. (|82|) - 
( |S5| ) is rather complicated and not very well investigated. 
There are several different approximations which can be 
used to solve these equations in the context of preheating. 
We will describe them in this section. 



X 




4>{t) 

b) 

FIG. 12. Two diagrams for the polarization operator of 
the field (p. Thin and thick lines represent the fields rj> and 
X respectively. The first diagram corresponds to the Hartree 
approximation which takes into account the contribution of 
(x 2 )- The contributions of these two diagrams to the effective 
mass of 0-particles can be comparable to each other. 

B. Hartree approximation 

The simplest way to take into account the backreaction 
of the amplified quantum fluctuations x is to use the 
Hartree approximation, 

4> + 3H<f) + m 2 <f) + g 2 (x 2 )<p = , (86) 
where the vacuum expectation value for x 2 is 

oo 

<X 2 ) = m / dkk2 \ X ^)\ 2 - (87) 
2-K z a 6 J 

o 

Quantum effects contribute to the effective mass of 
the inflaton field as follows: m| = m 2 + g 2 (x 2 )- The 
Hartree approximation corresponds to the first of the two 
diagrams of Fig. [l2[ 

Initially, we have no fluctuations fk{t), and we can 
use Eq. ( |38| ) for the modes One can express (x 2 ) in 
terms of the ak(t) and (3k{t) coefficients describing the 
resonance: 




24 



This formal expression may need to be renormalized. 
The WKB expansion of the solution of equations (0) 
provides a natural scheme of regularization p7[ . How- 
ever, in our case the coefficients ctk and (3k of the Bogoli- 
ubov transformation appear due to particle production 
(as opposed to vacuum polarization), so the integral in 
Eq. ( |S8| ) is finite and does not require further regulariza- 
tion. 

Let us estimate (x 2 ) from Eq. 
suits of the previous section. For 



ation of x particles we have \(3 k 



using the re 
the resonant ere- 



l„2fi k m t 
2 ' 



jdt 



|/3 fc | 2 cos(2 J ujdt - arga fc 



Re(a k p* k e 2l I L 

arg/3fe). For ui ~ 5</>(i) = sin mi the phase in this 
expression is equal to cos mt plus a small correction 
0(k 2 ). Due to this small correction, the term 2^ cos mt 
acquires a numerical factor C < 1 after the integration 
Jd 3 k: 



<x 2 ) 



1 + C cos 2g<I>cosmt 



27r 2 a 3 



dkk 2 



-rife . 



(89) 



In the broad resonance case when (f> > 4>* (i-e. for most 
of the time), one has - <C g<f>, to ~ g|(/>(t)|, and therefore, 



(x 2 > - ( 



I + C cos 



2gQ cos mt 



rn 



9W)\ 



(90) 



This means in particular that in the broad resonance 
regime the effective mass squared of the background field 
<f>(t) in the Hartree approximation 



1 + C cos I cos mt 

V m 



9n x 



(91) 



oscillates with two frequencies. One is the frequency of 
oscillation of |</>(t)|, which is equal to 2m. In addition, 
when (j>(t) ^ «3>, the effective mass squared m\ oscillates 
with a very high frequency ~ 2<7<I> 3> m. The amplitudes 
of both oscillations are as large as the maximal value of 
g 2 (x 2 }- One can easily identify both types of oscillations 
of (x 2 ) in the numerical simulations of Khlebnikov and 
Tkachev @. 

The resulting equation for the field 4>(t) looks as fol- 
lows: 



■ + 3H<j) + m 2 4> 
+ gn x (l + C cos 



2g$ cos mt \ (f> 
m )\$\ 







(92) 



The last term in this equation oscillates with a frequency 
~ 2<?<f>. In the broad resonance regime with g$ ^> m the 
high-frequency oscillation of this term does not much af- 
fect the evolution of the field <f>(t) because the overall sign 
of the term C cos 2 9 j '™ smt changes many times during 
each oscillation of the field </>. One may wonder, how- 
ever, whether these high-frequency oscillations may lead 
to a copious production of </>-particles. A preliminary 



investigation of this issue shows that the quasi-periodic 
change of the last term in Eq. (|92] ) does not lead to para- 
metric resonance, but a non-resonant particle production 
is possible because the effective mass changes in a very 
nonadiabatic way: ~ gm<§ ~ k 2 m 2 . 

In the first approximation one may neglect this effect 
and write Eq. (|92|) as follows: 



3H(f> 



m 2 (f> 



9 n x 



. 



(93) 



Even in this simplified form the last term of this 
equation looks rather unusual. It is not proportional 
to (f>, which would be the case if %-particles gave a <p- 
independent contribution to the effective mass of the field 
cf>. In our case this contribution is inversely proportional 
to |</)| . As a result, the field <j> behaves as if it were oscil- 
lating in the effective potential <?n x |0|. 

To estimate the change in the frequency of oscillations 
of the field due to the term gn x -^ in Eq. (|93|), one 

can neglect the term 3H(p in the equation for the ho- 
mogeneous field <fi, because 5 « m at the end of the 
first stage of preheating, when the term gn x jjj^ becomes 
important. Let us find when the frequency increase due 
to the interaction with ^-particles becomes greater than 
the initial frequency m. In order to do this one should 
solve the equation (f> = —gn x in the interval < <p < $. 
The time during which the field <j) falls down from $ to 



is At 



This time corresponds to one quarter 



of a period of an oscillation. This gives the following ex- 
pression for the frequency of oscillations of the field <p in 
the regime when it is much greater than its bare mass 
squared m 2 : 



W0 



2V2 



7710 W 7710 . 



(94) 



Here m^ is the value of the effective mass of the field cf> 
at the moment when <j>(t) = $ (the oscillations of (x 2 ) 
being ignored). Therefore to estimate the change of the 
frequency of oscillations of the scalar field 4> one can use 
the standard expression m 2 = m 2 + g 2 (x 2 ) for the ef- 
fective mass squared of the field 0, where by (x 2 ) one 
should understand its smallest value per period, which 
appears for 4>(t) = <&. This implies that the frequency of 
oscillations of the inflaton field does not change until the 
number of x-particles grows to 



?rr<I> 2m i 



g 



,1/2 



(95) 



This is a very important criterion which defines the du- 
ration of the first stage of preheating where the backre- 
action of the created particles can be neglected. 

For future reference we include here expressions for 
the energy density and pressure of the nonrelativistic x~ 
particles. The contribution of x-particles to the energy 
density p x {4>) °f the oscillating field 4> m terms of otkit) 
and Pk(t) is given by 
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where \j3k 



n k . This expression does not have any 



high-frequency modulations which we have found for the 
Hartree term ^ (x 2 ) 4> 2 ■ During the main part of each 
oscillation of the field <j), the field x has mass much greater 
than the range of the integration ~ k* , which means that 
w ~ and 



PxO) 



9[ 



2n 2 a 3 



dk k 2 nk = g\(f>\n x 



(97) 



The contribution of x-particles to pressure in terms of 
a k (t) and /?*.(*) is given by 



dkk lj 



Re (a t f)Je- 2 ' £"'") + 
2g$ cos mi 



-g\cj)\n x C cos 



m 



The last equality holds in the nonrclativistic limit, for 
S> 4>*. Average pressure in this regime is equal to zero, 
as it should be for nonrelativistic particles. 



C. Is the Hartree approximation sufficient for the 
calculation of particle masses? 

In the previous subsection we investigated the change 
of frequency of oscillations of the classical background 
field 4>(t) due to its interaction with x-particles, see Eqs. 
( |93| ) and (|94|). What about the spectra of perturbations 
8<jf! In order to answer this question one should cal- 
culate both diagrams shown in Fig. [l^. The first of 
these diagrams, Fig. |l2|a, gives the same contribution 
11^ = g 2 (x 2 ) as the one which we already calculated 
when we studied oscillations of the field <fi(t). As we 
have seen, in the situation where fluctuations Xk{t) are 
amplified by resonance, even the calculation of this sim- 
ple diagram is rather nontrivial and leads to an unusual 
result ( P2|). The calculation of the polarization operator 
n 2 , Fig. |l2[b, is much more involved. Similar diagrams 
have been ignored in all previous papers on preheating. 
Let us try to understand, however, whether II^ can be 
neglected as compared with 11^. A positive answer to this 
question would imply that the Hartree approximation is 
sufficient not only for the investigation of the oscillations 
of the field <f)(t), but also for finding the spectrum of per- 
turbations of the field <f). 

Usually when one calculates similar diagrams at high 
temperature, the polarization operator ilj, in the high- 
temperature limit is proportional to T 2 , whereas is 
less divergent at large momenta and therefore grows only 
as T. Therefore in the high-temperature approximation, 



the first diagram, which corresponds to the Hartree ap- 
proximation, gives the leading contribution. In our case 
this issue should be reconsidered because the leading con- 
tribution to the diagrams is given by particles with large 
occupation numbers and relatively small momenta. 

The backreaction of created particles becomes essen- 
tial only at later stages of reheating, when, as we will see 
shortly, H <C m. Therefore at that stage one can neglect 
the expansion of the universe when calculating polariza- 
tion operators, and it is more convenient to perform all 
calculations in terms of the usual, physical (rather than 
comoving) momenta k and the modes Xk(t). Therefore 
throughout the rest of the paper we will use physical 
momenta, k, p, etc. During the last stages of reheating 
they remain almost constant, but in order to relate them 
to the original physical momenta for each mode Xk one 
should remember that physical momenta are redshifted 
as a _1 (t). 

To calculate lii one needs to know the Green function 
of the field x i n & n external field (j>(t), which is given by 



G x (x,x')= / d 3 kT[ Xk (t)xk(t')l 



ik(x— > 



(99) 



where T stands for time-ordering. The calculation of the 
diagram for n 2 , Fig. [l2]b, using this Green function for 
the internal lines of the field x is rather tedious. There- 
fore, we will make certain simplifications. Consider the 
broad resonance regime q 3> 1 at a time when <p(t) 3> </>*. 
At this stage there is no particle production, and the adi- 
abatic form (p2[) can be used for the eigenfunction Xk(t). 



Consider a time interval At < 



near the time when 



the inflaton field (f>(t) reaches its maximum, $. During 
this short interval, one can neglect the expansion of the 
universe and the change of the field (j>(t), i.e. one may 
take <p(t) w <f>. The Green function in the space-time rep- 
resentation consists of two parts. The first part is similar 
to the standard Green function in Minkowski space in 
the fixed background field <p. The second part contains 
the high frequency modulation e lul ^ t+t \ Both terms are 
of the same order. One can show that in this regime the 
first term in the expression for the Green function ( |9Sj ) 
has a simple form in the momentum representation: 



G(k) 



k 2 - 



— + 2nn k S(k 2 - m 2 ,) 



(100) 



Here m x = g<fi(t), and k is a physical momentum. The 
first term in this equation is the standard Green function 
for quantum fluctuations in the vacuum. The second 
term is proportional to the occupation number n k = \Pk\ 2 
of the x-particles. 

The second part of the full Green function containing 
the modulation e tuJ ( t+t ) does not have a simple interpre- 
tation in the momentum representation. Omitting this 
part does not affect the order-of-magnitude estimate of 
the polarization operator. This can be most easily seen 
for the diagram Fig. |l2]a, where the calculations are much 
simpler. Indeed, with the complete Green function (|99|) 
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one can immediately reproduce the result (|9^) fo r th e 
diagram in Fig. [l2|a. Meanwhile, if one uses Eq. ( 100 ), 



then in the large n& limit one gets the first, nonoscillating 
term in the brackets of (|9l|): 



ni 



9 
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d p2~K5{p 2 — m?,)n p = 



(101) 



The part of the Green function containing the modula- 



in this case gives us the second (rapidly 



tion e luJ (*+*') 
oscillating) term in Eq. (91). 

Thus, whereas in the first approximation one can in- 
terpret the growing modes of the field \ during paramet- 
ric resonance as normal particles on the mass shell with 
the standard Green function (100), this interpretation in 
general is not quite adequate and may lead to the loss of 
some terms such as the oscillating term discussed above. 
Still we correctly reproduced the most important part of 
the polarization operator nj,. 

Let us try to estimate the polarization operator n| 
using the simple Green function (|l00| ) for \<j>(t)\ w $. The 
general structure of the polarization operator is given by 



. 9^ 
l ~2 



J d 4 pG(p)G(p-k±q) 



(102) 



The sign of q depends on whether the external field <j){t) 
brings the momentum go = m, q = to the two ver- 
tices of the polarization operator or takes this momentum 
away. 

It is not our purpose now to perform a complete calcu- 
lation of LT^ in this paper because we do not need to know 
the exact spectrum of perturbations 5<f>. Our main goal 
here is to find out whether or not LT^ may contain terms 
comparable to the Hartree operator H^. And indeed, if 
one calculates, for example, the diagram where the ex- 
ternal field <j)(t) brings a momentum qo — m, q = to 
the first vertex and takes it away from the second vertex, 
one finds (ignoring factors 0(1)) that this contribution 
to the real part of n 2 for fco = m, k = in the limit 
rip 3> 1 has the same structure as IT^: 



Reffi 



g 4 * 2 
(2tt) s 



i p d 3 p 



Po (Pa' 



^ . (103) 

q> 



Here po = uj = \/p 2 + g 2 $ 2 w g$> for a typical reso- 
nant mode with g 2 ® 2 > p 2 ~ gra$ > m 2 . Thus, for 
\4>{t)\ s» $ the second polarization operator of Fig. [l2| 
contains terms of the same order of magnitude as the 
value of the polarization operator in the Hartree approx- 
imation. This result indicates that one may need to go 
beyond the Hartree approximation used in many papers 
on preheating. 

This result looks paradoxical. In particular, one could 
argue that the Hartree approximation is closely related 
to the 1/N approximation, which is expected to give ex- 
act results in the limit N — > oo. Indeed, instead of a 



single x-field one can take N fields Xi with the interac- 
tion ^y0 2 xf- The Hartree diagram is proportional to 
g 2 , i.e. it survives in the limit N — » oo , whereas the ex- 
pression for the polarization operator is proportional 

4 

to jj. That is why usually at large N one can neglect 
contributions like n 2 as compared with n^. Indeed, this 
would be true in our case as well if the field x had a 
large ^-independent mass. But in the theory we are dis- 
cussing now its mass squared is jj(f> 2 - As we have seen, 
when one calculates Tii this mass squared appears in the 

4 

denominator. As a result, the factor jr in front of the 
diagram becomes g 2 , so that this diagram also survives 
in the limit N — *■ oo and has the same order of magnitude 
as the Hartree diagram in the 1/N approximation. This 
means, in particular, that without a complete calculation 
of n 2 one cannot be sure that the 1/N approximation 
gives a correct spectrum of particles in the limit N — > oo 
when applied to the theory of preheating. 

To avoid misunderstandings we should reiterate that 
this problem appears in the calculations of the effective 
masses of the ^-particles but not in the calculation of 
corrections to the equation of motion of th e bac kground 
field (j)(t), which was our main goal in Sec. VIII. 



D. Classical approximation to the self-consistent 
dynamics 

Fluctuations of bose fields generated from vacuum by 
an external field in the large occupation number limit can 
be considered as classical waves with gaussian statistics, 
see e.g. jy|. Therefore in the first approximation all fields 
X, S(f> can be treated as interacting classical waves. This 
makes it possible to study preheating by investigating 
a system of nonlinear classical equations or by lattice 
numerical simulations of the interacting classical scalar 
fields HIlS 

The Fourier decomposition of the Klein-Gordon equa- 
tions of the interacting fields can be reduced to mode 



equations. The mode equation for Xk 



a 3/2 Xk is 



X k +\-^+g'r[t)\X 



(27r) 3 a 3 / 2 



d*k'X k _ k ,<p k , 



ir 1 d 3 k'd 3 k"X k _ k , +k „ip kl ip k „ 



(27ra) 3 

The mode equation for 5<f>k{t) = a~ 3 ' 2 (pk(t) is 
k 2 



(104) 
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(105) 
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The first term in the r.h.s. of this equation describes 
rescattering of x-P ar ticles on the classical field 4>{t), 
which leads to 0-particle production. The second term 
describes scattering of ^-particles and x-p&rticles. Cor- 
rections to the effective mass of the modes 4>k appear as a 
result of the iterative solution of the system of equations 
which we now present. 

The equation for the oscillating background field <j>(t) 
looks as follows: 



cj) + 3H(j) + m z 4> 

,2 



9 2 <t> 



d 3 k'X£, 
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The first term on the r.h.s. of this equation is propor- 
tional to the polarization operator 11^, which is shown in 
Fig. 12 a. The second term describes rescattering, which 



is related to the imaginary part of the polarization oper- 
ator II?, Fig. |l2|b. Neglecting this term, one reproduces 
Eq. ( pq ) with the term containing J d 3 k\X k \ 2 playing the 
role of the induced mass. Thus the classical approxima- 
tion reproduces the Hartree approximation, but it also 
takes into account effects related to rescattering. 

In the beginning one can neg lect <Pk{t) and the corre- 
sponding integral terms in Eq. (|l04j) . Later, the fluctua- 
tions Xk(t) are amplified by the resonance and giv e ris e 
to (pk(t) fluctuations via the integral terms in Eq. ( 105 ). 
When the amplitude of fluctuations <fk(t) grows signifi- 
cantly, they begin to contribute to the integral terms of 
Eq. ( |104p . We will show (see Sec. |X B[ ) that the ampli- 
tude lpk(t) grows with time as e 2,lmt . Therefore the num- 
ber of particles corresponding to 8(f) fluctuations grows as 
much faster than n v . The interaction terms in 



Eqs. (104) and (105) can be interpreted as scattering of 
X particles on the inflaton field. Because of the very fast 
generation of 5<f> fluctuations, \S<f)\ 2 oc e 4 ^ 4 , the process of 
rescattering can be very important. However, it is not so 
easy to evaluate its full significance for the efficiency of 
the resonance. For example, if the particles <p produced 
during rescattering have small momenta fc, they cannot 
be distinguished from the homogeneous oscillating scalar 
field, and therefore they do not make any difference to 
the development of the resonance, see the discussion of 
this issue in Sec. H Therefore we need to know not only 
how many 5<j) particles are produced, but also whether 
they are "hard" particles with large momenta or "soft" 
particles with small momenta. We will return to this 
question in Sec. KL 



IX. TWO STAGES OF PREHEATING, 
RESCATTERING BEING NEGLECTED 



Previously, we were mainly following the evolution of 
each particular mode \k- Now we will study their integral 
effect in an expanding universe. 



As we have found in the previous section, the develop- 
ment of broad parametric resonance can be divided into 
two stages. In the first stage n x < backreaction of 

the particles \ can be neglected, and the frequency of os- 
cillations of the field </> is determined by its mass m. (We 
will argue later that at this stage their scattering also 
does not lead to any important effects.) In the second 
stage n x > , and the frequency of oscillations of the 
field <f> becomes determined not by its bare mass, but by 
its interaction with x-particles. Now we will study the 
first and second stage of broad parametric resonance. 

We begin with the first stage when the backreaction 
of created particles can be neglected. Then we consider 
the second stage where backreaction is important assum- 
ing a certain hierarchy of the feedback effects: effective 
mass of the inflaton is changed first, and rescattering 
may become important afterwards. In this section we 
will neglect rescattering. In the next section we will dis- 
cuss rescattering and the validity of the assumption men- 
tioned above. 



A. The first stage of preheating: no backreaction 
and no rescattering 

In the first stage of preheating one can ignore the back- 
reaction of created particles on the frequency of oscilla- 
tions of the field <j>{t). As we have found in Sec 
this stage ends at the moment t\ when 
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In the next section we will show that the effects related 
to rescattering also do not alter the development of the 
resonance during this stage. In this section we will es- 
timate the duration of the first stage t%, the number of 
inflaton oscillations Ni at the time t = t\, the number of 
created particles n x {t\), the energy density of these parti- 
cles p x (ti) and the value of (% 2 (ii)). We will use symbols 
$, q and fc* without any indices for the running (time- 
dependent) values of the amplitude of the field (j>(t), of 
the (/-factor, and of yjgm<&{t), whereas, for example, q a 
will correspond to the value of q at the beginning of pre- 
heating, and qi will correspond to its value in the end of 
the first stage of preheating. 

One can use Eq. (|64|) to estimate n x . First one should 
determine which fluctuations Xfc are amplified during the 
entire period of the resonance. The fluctuations am- 
plified by the broad resonance have physical momenta 
k < fc*/2 ~ see Eq. (|I]). (More precisely, one 

may expect k < k^/y/n, see Eq. (pq).) Then the am- 
plitude <E> in this expression decreases as about M p /3mt. 
Therefore, the resonance width decreases as k ~ t^ 1 ! 2 1 
whereas redshift of the momenta of previously produced 
particles occurs as a -1 ~ <~ 2 / 3 , i.e. somewhat faster. 
(In terms of comoving momenta fc, the resonance width 
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grows as k ~ a(t)yJgm<&/2 oc i 1 / 6 .) This means that 
those modes which have been amplified at the first stages 
of the process continue to be amplified later on. There 
are modes which were outside of the resonance band in 
the very beginning, but entered the resonance band later. 
However, after a time ~ (2/im) -1 the fluctuations which 
have been amplified from the very beginning will be ex- 
ponentially larger than the "newcomers". Therefore the 
modes which do not enter the resonance band from the 
beginning typically give a subdominant contribution to 
the net effect. 

Thus, with reasonably good accuracy, during the first 
stage of preheating one may consider only those fluctu- 
ations which have been amplified from the very begin- 
ning. This is important because it means that in all 
integrals one should consider only momenta which ini- 
tially, when a(to) = 1, Q(t) — $o, were in the interval 



fc(*o) < fc* (*o)/2 ^ ~ m q f /V2. (108) 

where on = ^r%- 

The most important element of our calculations is the 
exponentially growing occupation number of particles 
with k < k m : n(t) cx e 2fimt . Here \x is an effective in- 
dex which describes an averag e rate of growth for modes 
with k < k*, see Sec. VII C . In our model \x depends 
on g, but not very strongly, see the table in Sec. jv|. 
Typically it is in the range 0.1 — 0.2. For definiteness, in 
our estimates we will use fi — 0.13 which we have found 
numerically for a certain range of values of the coupling 
constant g, see the table in Sec. VI. As we will see, in 



the context of our approach an error in our estimate of /x, 
say of 10%, does not create an exponentially large error 
in the final result (contrary to the remark of jL7|); it only 
leads to an error of 10% in the calculation of the dura- 
tion of the first stage of preheating. Our final results will 
be even less sensitive to the value of the subexponential 
factor in Eq. @. 

Substituting Eq. (|l0g|) into Eq. <Mf, we find 



<(*) 



(smf ) 3/2 



O 2jum(t-io) 



647r 2 a 3 \J ir/j,m(t — to) 
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where to is the beginning of the inflaton oscillation. The 
convention we used in Sec. VI is that <o — 7t/2to, which 



gives $o — Mp/5 and qo — 10 g . Our choice is also 
very close to the convention of ref . . (This particular 
choice is not going to be important because the total 
duration of the process is much greater than to.) With 

'2mt\ 2 / 3 



this choice of to we have a(t) 
has| 



' . For<>t one 



n x (t) ~ 10 



_ 4 (gmM p ) 3/2 
(mt) 5 /2 M i/2 



(110) 



Now we have to sub stitute Eq. ( 110 ) and $(f) ~ 
M p /3mt into Eq. (107). The result can be transformed 
into an equation for t\\ 



ti 



1 , 10 6 m(mii) 3 
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An approximate solution of Eq. (Ill) for \i ~ 0.13 
is t\ ~ In gijj" @- As we will see soon, this is a 
good estimate not only for the duration of the first stage 
of preheating, but for the duration of the whole process, 
because the second stage of preheating typically is rather 
short. 

For a realistic value m ~ 10 _6 M p in chaotic inflation 
in the theory m 2 </> 2 /2, our estimate gives Q 



5 , 15 
In — 

4/xm g 



(112) 



For instance, for fj, = 0.13 and g = 0.1 one has t\ ~ 
50m -1 ; for g = 10~ 2 one has t\ ~ 70m _1 ; for g = 10 -3 
one has t\ ~ 90to _1 , etc. 

The value of the field $i = at the end of this 

first stage is given by 



M p iuM v , 10 12 m 
— = hi — e ■ 

3mti 3 g M p 



(113) 



Another important quantity is the value of the parameter 



q = at the end of the first stage: 



1/2 3*1 2g^M p _j 10 12 m 
ffi = = in — 

yi 2m 3m g 5 M p 



(114) 



*The total duration of the first stage of preheating in our 
model typically is about 10 2 m _1 . If one compares the red- 
shifted value k*(to)/a(t) of the physical momentum which was 
equal to fe*(to) ~ \/gm$o at the beginning of preheating, and 
the running value of fc»(t) ~ y/ gm$(t), one finds out that 
in the beginning these two quantities coincide, whereas at 
the end of the first stage of preheating the running value of 
fe* is greater than the redshifted one by only a factor ~ 2. 
Therefore at the end of the first stage and at the beginning 
of the second stage of preheating instead of calculating the 
redshifted value of k*(to)/2 one may simply use the condition 
fc < fe*(t)/4 for the growing modes. 



^Equation (110) is a starting point for our further esti- 
mates. To derive this equation we used the theory of suc- 
cessive para bolic scatterings. However, the general structure 
of Eq. (110) can be easily understood even without any use of 
this theory. As we already mentioned, the value of can be 
obtained by solving the Mathieu equation numerically in an 
expanding universe, see Sec. VI. One can make a simple esti- 
mate of n even without using a computer. Indeed, we know 
that the parameter along the line A — 2q changes from 
to 0.28 jj]. An average of these two numbers, 0.14, provides 
an excellent approximation to the true value of p. 



29 



To find the typical occupation numbers at the end 
of the first stage of reheating, let us remember that 
n x = / k 2 dknk, and that integration typically goes 
from to the physical momentum ~ fc*/2. This gives an 
estimate 



487r 2 n 



v 



(115) 



The occupation numbers of x-particles rife(ti) by the 
end of that stage can be estimated as nkitx) — 3 x 
10 W 1/4 , see Eqs. (§TJ) and (jTTTJ). 

Using the results of this section, for different values 
of the coupling constant g one can estimate the initial 
value qo of the parameter q, its value q\ at the end of the 
first stage of preheating, the value $i, and the number of 
oscillations N± which the field makes from the end of 
inflation to the end of the first stage. In the table below 
we give somewhat rounded numbers: 



g 


9o 


9i 




Nx 


10~ a 


10 4 


3 


3.5 x 10" a 


15 


io- y 
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550 


5 x 


11 
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10 8 


10 b 


7 x 10~ J 


8 



by 



The energy density at the end of the first stage is given 



! 2 $2 8/z 2 m 2 M 2 _ 2 10 12 m 
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It is worth comparing the frequency of the inflaton oscil- 
lations to with the Hubble parameter at that time: 

/F2$i 8u _! 10 12 m 

For instance, for /i = 0.13, g — 10~ 2 , to = 10 _6 M P one 
has 



J?(ti) - 10" 2 to. 



(118) 



Thus, at the last stages of preheating (though not at the 
beginning) one can, in the first approximation, neglect 
the expansion of the universe. 

At that time, when g 2 (x 2 ) — m2 , the total energy 
density (on the r.h.s. of Eq. (|82|)) becomes approxi- 
mately equally distributed between the interaction en- 
ergy V x {4>) — g^\n x = m 2 <&\ and the potential energy 
density to 2 $ 2 /2 of the field cj). The kinetic energy of 
X-particles can be estimated as ((Vx) 2 ) — k 2 (x 2 ) — 
5$ito(x 2 ) ^ m 2 $ 2 ^ ~ to 2 $ 2 g~ 1/2 . 

If preheating does not end with the end of the first 
stage, i.e. if q\ >• 1/4, then the kinetic energy remains 

small: ((Vx) 2 ) — to 2 * 2 g^ 1/2 < g 2 $\{x 2 ) — m 2 <&\. 
However, if at the end of the first stage q% ~ 1, then at 



that time a considerable fraction of the energy of the in- 
flaton field will have been transformed into the kinetic en- 
ergy of the x-particles: ((Vx) 2 ) — to 2 $ 2 q x 1 ^ 2 ~ to 2 $ 2 . 

Let us find the range of values of the coupling constant 
g for which preheating ends during the first stage and 
for which investigation of backreaction is not necessary. 
Without taking account of the backreaction preheating 
ends at the time tf when g&(t) drops down to to, which 
gives tf ~ (see Sec. VI). Therefore, preheating ends 



in the first stage if tf ^ t\, i.e. if 



/ 

. 4m , 15 
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For our values of parameters this gives the condition g < 
3 x 10 -4 . In our convention, this corresponds to an initial 
value q Q < 10 3 . 

In this regime the total number density of x-particles 
created during preheating is given by 



gM p 



exp • 



2g[iMp 
3 m 



(120) 



and the x-fluctuations at the end of this stage are given 
by 



(x 2 



gMp 



exp ■ 



2gjjMp 
3to 



(121) 



Eq. (|119| ) implies that for g rss 3 x 10 4 this quantity 
should coincide with the value of (x 2 ) at the end of the 
first stage of preheating, (x 2 ) = r -^r- Thus, for g sa 3 x 
10 -4 one has 



3 x 10 ib GeV 



(122) 



The possibility to obtain enormously large fluctuations 
of the field x is one of the most remarkable features of 
preheating. For comparison, if the field x were in a state 
of thermal equilibrium, the dispersion of its fluctuations 
would be given 

by = T/2V3 @. Therefore in 

order to obtain yfl^ 2 ) - 3 x 10 16 GeV one would need 
to have T > 10 17 GeV, which is practically impossible 
in the context of inflationary cosmology. Here such fluc- 
tuations can be generated prior to thermalization due to 
t he r esonance at the stage of preheating. Fluctuations 
(122) change the effective masses of particles interacting 
with the field x- The simplest way to study this possi- 
bility is to add to our model another scalar field r\ with 
a potential describing symmetry breaking, for example, 
ViViX) — ^K 7 ? 2 — r /o) 2 + 7 7 2 X 2 ]- F° r sufficiently small 
A this addition does not affect preheating and does not 
change any of our results concerning (x 2 ) • It is obvious 
that the generation of perturbations (x 2 ) leads to sym- 
metry restoration in this model for r/o \J (x 2 ) 011 a 
scale up to 770 - 10 16 GeV |JTJ. Such effects may have 
important cosmological implications 43 1 . 

Thus, we can distinguish between different scenarios 
depending on the coupling constant g. 
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For g <C 3 x 10~ 4 the broad resonance ends during the 
first stage. In this case parametric resonance is not effi- 
cient enough to transfer a significant part of the energy of 
the inflaton field to the energy of x-particles. The most 
important part of the process of preheating in such the- 
ories is described by the elementary theory of reheating 



For g ~ 3 x 10 , at the end of the first stage 
<7i <~ 1/4, and the energy becomes approximately equally 
distributed between the energy of the oscillating scalar 
field <p & n d the energy of ^-particles produced by its os- 
cillations. 

For g > 3 x 10~ 4 the broad resonance continues after 
the end of the first stage. To investigate the further de- 
velopment of the resonance one should study quantum 
effects which could be produced by the ^-fluctuations in- 
teracting with the oscillating field 4>(i). 

Before doing so, let us remember that the presence of 
the interaction g 2 4> 2 x 2 typically leads to radiative cor- 



:',2tv- 



rections to the effective potential of the type 
For g > 10 

2M- for 6 



ln< 



" 3 this term becomes greater than the term 

~ 4M p , when the density perturbations de- 
termining the structure of our part of the universe were 
produced. Thus one may argue that in models of the 
type considered above g should be smaller than 10 -3 . If 
g belongs to the narrow interval between 3 x 10~ 4 and 
10~ 3 , reheating ends soon after the end of the first stage, 
and the effects of backreaction are only marginally impor- 
tant. For the description of preheating in such theories it 
is sufficient to use the simple estimates obtained in this 
section. 

However, in supersymmetric theories radiative correc- 
tions from bosons and fermions have a tendency to cancel 
each other. In such theories the coupling constant g can 
be much greater than 10 -3 . Therefore, we will continue 
to consider all possible values of the coupling constant g 
without assuming that g < 10~ 3 . 



B. The second stage of preheating, neglecting 
rescattering 

We defined the second stage of preheating, t > t%, as 
the stage when the frequency of inflaton oscillations due 
to the feedback of amplified x fluctuations is no longer 
m but is determined by the backreaction of x-particles. 
In this section we will investigate preheating neglecting 
rescattering. In the next section we will discuss the valid- 
ity of this assumption. The frequency of the infla ton os- 
cillations during this stage was derived in Section VIII B 
and given by Eq. (p4[). Since this frequency is much 
greater than the bare mass m, the second stage is much 
shorter than the first one. Indeed, at this stage each os- 
cillation takes a time which is much shorter than 2itm~ 1 , 



whereas the number of particles, as before, grows as 
e i7T ^ N , where N is the number of oscillations. Therefore 
the number of particles can grow exponentially within a 



time which is much shorter than H~ l . This implies that 
one can neglect expansion of the universe and the corre- 
sponding decrease of the total energy density of matter 
during the second stage of preheating. 

Let us consider the inflaton oscillations 6(t) during the 
second stage. From Eq. (|93| ) in the limit H <C we 
have 



(j) + m 2 <j) + gn x sgncj) = 



(123) 



where sgn<^ is ±1 depending on the sign of the value </>, 
n x (t) is a function of time, the expansion of the uni- 
verse is neglected, and m 2 4> <C g 2 n x sga(f>. The solution 
of this equation 6(t) consists of a sequence of segments of 
parabolas with opposite orientation that are symmetric 
relative to the t-axis and match at 6 = 0. The equa- 
tion for the modes Xk(t) will contain the square of g<j){t) 
instead of g 2 <& 2 sin 2 mi. Thu s, the behavior of Xk{t) f° r 
4>{t) determined by Eq. (122) is somewhat different from 
the behavior of Xk as described by the Mathieu equation. 
Nevertheless, this is not a real problem here. 

Indeed, if one does not take backreaction into acco unt, 
then, according to our investigation in Sec. VI, the 



system spends half of the time in the broad resonance 
regime, and another half of the time in the regime with 
q ~ 1, so this regime is very important. However, let 
us consider the effects of backreaction. The parameter 
q = g 2 ® 2 /4m 2 at the second stage can be estimated us- 
ing the "effective mass" (or, more exactly, the frequency 
of oscillations of the field <j>) m 2 ^ ~ gn x /§ ~ g 2 (x 2 ) @)- 
This gives q ~ g$ 3 /4n x . The end of the resonance, as 
before, occurs at q ~ 1/4, see below. The number of x 
particles grows exponentially, so during the previous os- 
cillation one had q ~ e 47rM /4 ~ 1, and during the previous 
oscillation q was much greater than 1. Therefore during 
all the time except the last one or two oscillations the pa- 
rameter q was very large, the resonance was very broad, 
and it could be described by the theory of stochastic res- 
onance. This theory is very robust; it depends only on 
the speed of the field (j> near 6 = 0. Thus, the difference 
between the Mathieu equation and the equation for the 
modes Xk m the field 4>(t) satisfying Eq. ( |123] ) in this 
context becomes unimportant. 

On the other hand, at the time when q decreases, the 
structure of the first resonance band becomes important. 
We investigated this issue by solving equations for the 
modes Xk n umer ically. We found that if the field 6{t) 
obeys Eq. (123), the structure of the first resonance 
band for Xk at small momenta is very similar to that 
of the Mathieu equation. Therefore, the second stage of 
preheating in this case ends when 
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1/4, 



(124) 



just as before. This happens at some moment ty. when 
5 $ 2 ~ m^h), $ 2 ~v/0t% n x (t 2 )~^^. (125) 
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At this time the total energy density becomes approxi- 
mately equally distributed between the kinetic energy of 



X-particles 



(x 2 ) and the energy ~ g&n x of their 



interaction with the field (p (which includes the potential 
energy of the field <j>). This energy should be equal to the 
total energy of the system at the time t\ , which is given 

by 3m 2 $1 . The final value of the inflaton field at the end 
of resonance is 



$ 2 



1/4 



(126) 



Thus, $2 is somewhat smaller than $i for qi > 1: 
2 V 3. 9 g 5 M p J 



(127) 



To find the typical occupation numbers of the modes 
with k ~ fc* at the end of the second stage of reheating, 

let us remember that — T 3 " x ■ This corresponds to 
enormously large occupation numbers pj 



n k (t 2 ) ~ 10 2 5 



2 „— 2 



(128) 



This result indicates potential problems with the per- 
turbative investigation of preheating at the end of its sec- 
ond stage. Adding extra internal lines of the diagrams 
may introduce enormous factors ~ 10 2 g~ 2 , which 
may cancel extra degrees of g 2 which appear in the higher 
order corrections. 

In order to calculate the duration of the second stage 
let us note that n x {t 2 ) ~ n x (ti)e 47r/xAr2 . One can show 



that 



™x(*i) 

ond stage is 



1/4 



Therefore, the duration of the sec- 



In 4qi 



1/4 



(129) 



Using the table of values of q\ given in the previous sub- 
section, one can conclude that the second stage may take 
from 2 oscillations (for g = 10~ 3 ) to about 10 oscilla- 
tions (for g = 1CP 1 ). Because of the growth of the effec- 
tive mass of the inflaton field, each oscillation takes much 



smaller time than — , so Eq. (Ill) for the duration of 
the first stage of preheating gives a good estimate for the 
total duration of the stage of broad resonance Q . 

Numerical estimates of $2 ~ \J (x 2 ) 2 show that it can 
be in the range of 10 15 to 10 16 GeV. As an example, 
for g = 1CP 2 , which corresponds to qo ~ 10 6 , one has 
i 10 16 GeV. An interesting feature of Eq. 



(127) is the inverse dependence of yj (x 2 ) 2 on ^he value 



of the coupling constant. 

Note that in addition to the high-fre quency oscillations 
with frequency ~ 5$ d iscus sed in Sec. |VIII B| , the ampli- 
tude of fluctuations yj (x 2 ) experiences oscillations with 
a frequency 2m. At the end of the second stage these 
two frequencies coincide. In all our estimates we cal- 
culated the minimal value of yj (x 2 ) which occurs when 



\4>(t)\ — It was important for us because this is the 
time which determines the frequency of oscillations of the 
field (pit). Near (pit) = the amplitude of fluctuations 
J (x 2 ) is greater than at \4>{t)\ — but close to the end 
of the second stage of preheating this difference becomes 
less significant. 

The results of numerical calculations of yj (x 2 ) per- 
formed in pof are in agreement with our estimates for 
g < 3 X 10~ 4 but give a few times greater value of yj (x 2 ) 2 
for larger g. The difference can be interpreted as a result 
of rescattering of x-particles during the second stage of 
preheating. 
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X. RESCATTERING 

Theoretical considerations contained in 
numerical lattice simulations of preheating 
dicate that there is another effect which should be incor- 
porated into the preheating scenario. In the context of 
the model investigated in this paper, one should consider 
the generation of inflaton fluctuations 8<p due to the in- 
teraction of x particles with the oscillating inflaton field 
(pit), and subsequent interaction betwee n x an d Sep fluc- 
tuations. We already discussed in Sec. VIII D the pos- 



sibility to describe this process by equations for classical 
waves. One may also represent the classical scalar field 
as a condensate of ^-particles with zero momentum, and 
interpret ^-particle production as a result of rescatter- 
ing of Y-particles and the ^-particles in the condensate 
[p4| , p5 29 3(]]. This "particle- like" interpretation of the 
interaction allows one to use the concept of cross-section 
of the interacting particles, and the Boltzmann equation 
for the occupation numbers. 

The theory of this process is rather complicated, and 
its interpretation in terms of the rescattering of ele mcn- 
tary particles is not universally valid, see Sec. [XB| . Still 
we can formulate the following apparently general results. 
First, there is a significant generation of rapidly growing 



fluctuations 



» 2 '"™* t due to the interaction between 



X-particles and the oscillating field <p(t). Second, the gen- 
eration of large fluctuations of 8<p can terminate the reso- 
nant creation of x particles only at the end of the second 
stage of reheating. In this section we will try to justify 
these statements. 



A. Generation of ^-particles by rescattering 

To evaluate the effects of rescattering we will use here 
the "particle-like" interpretation of rescattering. First, 
one should make an estimate of the cross-section a for 
the scattering of x particles with an effective mass g<p(t) 
and a typical physical momentum ~ fe*/2 = y/gmQ/2 on 
(j) particles of mass m with zero initial momentum which 
constitute the oscillating field (pit). The effective mass of 
the field x is time-dependent. This makes investigation 



32 



of their scattering rather complicated. However, in the 
broad resonance regime during the main part of the oscil- 
lation (for |0| > ~ i<&g~ 1/4 (fH)), the field x changes 
adiabatically. During this time, the effective mass of the 
field x &l so changes adiabatically, so one may consider 
%-particles as ordinary particles with an effective mass 
gcf>(t). We will also consider the oscillating scalar field 
4>{t) as a collection of particles c/> with an effective mass 
nitf, and number density = m^ 2 /2. 

We will suppose now that in such situation one can use 
the standard result for the cross-section for elementary 
particles (j) an d X with constant masses: 



(—) 

\dilJ CM 



teTpE4,E x (E4, + E x )\v 4 - v x 



(130) 



Here all energies , E x and velocities , v x are given 
in the center-of-mass (CM) frame and refer to the initial 
state, except for which refers to the final state. A4 2 
is the square of the matrix element, which is given by g 4 

During most of an oscillation one has \(j)\ > <p* ~ 
I$g- 1 / 4 ) an( j mx — gcj) ^ k* ~ y/gm^ffr. In this case 
both the ^-particles and x-particles are nonrelativistic. If 
one goes to the CM frame one finds that the <fi particles 



have a small speed 



^> v Y . Thus Ea 



2\/ 9^ 

j x ~ gcp. For g(f> 3> m$ the absolute value of the momen- 
tum of the ^-particles does not change after scattering, 

\p^\ k, ^y-^J 1 -^ <C m^. This gives, after the integra- 
tion of Eq. ( |130[ ) over <if2, a single particle cross-section 

4 2 

Now one should take into account that the actual cross- 
section will be much greater because the scattering occurs 
not in a vacuum, but in a state which already contains 
many bosons <f> and x- There are many x-particles from 
the resonance and many inflaton particles (p. Naively one 
would expect that the cross-section should be propor- 
tional to the product of the occupation numbers and 
n* in the final state. However, the corresponding terms 
disappear in the collision integral in the Boltzmann equa- 
tion, which takes into account all the channels of scatter- 
ing. Therefore in the investigation of enhancement of 
the cross-section due to the large occupation numbers 
of particles in the final state, one should consider terms 
proportional either to or nf, but not to n^n^. In the 
beginning of the process n x ^> n$, and the cross-section 

(7i should be multiplied by ~ 487r fc "*^ ■ This gives, 
for fa) 

Using this result, one can estimate the time for each x~ 
particle to experience one scattering with a 0-particle be- 



longing to the oscillating field <f>(t): r 

In particular, at the end of the first stage, n x 
which yields 



-1 1/2 

m 4> 1i 



0.5*-. 

n x 
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For g ~ 10 -3 this time is of the same order as the 
time of one oscillation of the field </>, see the table in 
Sec. IX. However, just one oscillation before the end of 



the first stage the density of particles was much smaller 
and rescattering was inefficient. For g > 10 -2 this time 
is much greater than the time of one oscillation, which 
means that rescattering occurs only during the second 
stage of preheating. 

In the "particle-like" picture the number of x particles 
does not change in each act of interaction (apart from 
its growth due to the resonance), but each collision re- 
leases one ^-particle from the homogeneously oscillating 
field <fi(t) Since the scattering time for each x-particle 
t oc n x (t), one may conclude that the number of free </>- 
particles grows with time as cx 5n 2 ,/$ 2 m0 cx e 4 *""* 1 . 
However, the true dependence is more complicated be- 
cause during each interaction the x~P ar ticles will slow 
down. This affects their subsequent interactions. 



B. On the validity of the "particle-like" 
interpretation of rescattering 

In the previous subsection we considered rescattering 
of particles during time intervals when <f>(t) > 4>*. At that 
stage x particles are nonrelativistic. In contrast, during 
the short time intervals At* ~ fc" 1 , when \(j){t)\ < (f)*, %- 
particles are ultrarelativistic, and their effective mass gcf> 
is very small comparing to their typical momenta ~ fc*/2. 
If one uncritically repeats the calculation of the rescatter- 
ing for the case of ultrarelativistic x particles in the time 
interval At* , one obtains a much higher cross-section and 

than that 



3?r n> 



a much shorter rescattering time t 
of the non-relativistic case of the previous subsection. 

However, within the very short time interval At* ~ 
g _1 / 4 m7 , one cannot use the standard methods of cal- 
culation |}2j developed for the investigation of processes 
which begin at t = — oo and end at t — +oo. The uncer- 
tainty principle tells us that during the time At* one can- 
not specify the energy of particles with an accuracy bet- 
ter than fc*. Therefore during the short interval At* one 
cannot tell the difference between a ^-particle with mo- 
mentum k = 0, belonging to the classical field 0(t), and 
a free ^-particle with momentum k < fc*, i.e. one can- 
not tell whether scattering occurred or not. This ques- 
tion can be answered only by observing the system for 
a longer time, comparable to rn^ 1 , but during the main 
part of such intervals the effective mass of each ^-particle 
is large, and cross-section is much smaller than the cross- 
section which one would obtain by naive application of 
the S-matrix approach during a small interval At*. In 
other words, we cannot use the standard formalism of 
particle scattering to describe scattering around zeros of 
the inflaton field. Another element missing in this formal- 
ism is that the field x is n °t m an n x -particle quantum 
state, but is in the squeezed state. (We have discussed 
already one of the nontrivial consequences of this fact, 
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namely the high-frequency modulation of (x 2 )-) Thus 
one may wonder whether one can trust the results of our 
calculations for the more safe situation when (p > , and 
what we can say about the contribution of the intervals 
with <p < 4>* to the net rescattering effect? 

Here we will outline a possible way to answer this ques- 
tion. Let u s co nsider the self-consistent set of equations 
( |104| ) and (105) for t he in teracting fields in the classical 
approximation. Eq. (105) describes the evolution of the 
S(j)k(t) fluctuation. Let us concentrate on the first inte- 
gral term in Eq. (105), assuming for the moment that the 
second term is subdominant until S(f>k(t) increases suffi- 
ciently. What we obtain is the equation for the forced 
oscillations of S(f>k(t). The homogeneous part of this inho- 
mogeneous linear differential equation has a simple Green 
function oc sin 17, k (t — t ' ) , where fl 2 . = k 2 + to? . Then the 
solution of Eq. (|l05|) with only the first integral term is 



<r 



dt' sinfi k (i - t')<f>(t') 



(2^)3fi k JQ 
cPk'x^{t') X v{t') + h.c. . (132) 



Here, as before, A: is a physical momentum. This solution 
expresses the function 5(j>k(t) via the known functions 
(f)(t) describing the inflaton os cilla tions, s ee E q. (||), and 
the functions Xk{t), see Sec. VTL Eq. ( |l32| ) takes into 
account all the regimes of 4>(t), as well as the resonant 
amplification of Xk- In particular, from this it follows 
that the amplitude (fik(t) grows with time as e 2Mm **, be- 
cause the amplitude Xk grows as e'™'**. Therefore the 
number of particles corresponding to 8<j) fluctuations is 
proportional to n 2 (i) oc e 4 ^™** , i.e. grows much faster 
than n Y . Another specific prediction which follows from 



Eq. (132) is that the random field (5</>(£,x) at the early 
stages of its generation will have non-gaussian statistics 
in contrast to the random gaussian field Xt.x - 

Let us further investigate the solution ( J132[ ) . The inner 
integral J d 3 k'xk-k'(t')xk'(t') is time-dependent. It is 
convenient to change the order of integr ation of J dt' and 
J d 3 k' . Then the r.h.s. of the solution (132) will contain 
terms like 



g 2 $e 



d 3 k' 



dt'- 



x e 



o V^k- k ' (t')u} k > (f) 

iSlki'+im^i' — i J dt" (t")+i J dt"uj k i(t") (133) 



where uJ k (t) = k 2 + g 2 <p 2 (t). During each half of the 
oscillation /3 k is constant, see Section [VII A . (Note that 
the coefficients /3 k corresponding to the classical waves 
will be dimensionless if one uses discrete modes k.) 

It is easy to see that different choices of signs in Eq. 
(133) correspond to different channe ls of scattering be- 
tween 5<j), X and <f> waves. The terms (133) correspond to 
the generation of 5(f> fluctuations due to the scattering of 
X and 4> waves. Obviously, one can leave in the inner in- 
tegral J Q only the segment J t (where t <tj + ^-) from 



the most recent cycle of the inflaton oscillation, when 
/3 k / is the largest. During this interval /3 k / is constant. 
Therefore to further investigate the inner integral f.* , we 
shall consider the variation of the phase of the exponent 
in Eq. (|l3|) 9 ~ -fl k t' + m^t' - f dt"w k - k >(t") + 



f dt" ui k i(t") within this time interval t — tj < —, 
which describes the interference of the four interacting 
waves <f>(t), <5</> k , Xk' an d Xk- k '- Earlier we estimated 
the integral /*' dt"u k ,(t") w ^ cos m i' + 0(k 2 ), see 
([77]). The crucial observation is that for the process 
Xk'<fto —> <50 k Xk'- k the large terms cosm^t in the 
expression for 9 are cancelled and the phase 9 does not 
oscillate within each half of the period, t — tj < —. 

As a result, the integral / dt cannot be reduced to the 
usual delta-function S (— f2 k + — w k _ k ' + u> k >), as one 
would expect in the "particle-like" picture. Instead, in 
the wave picture we will have nonvanishing contributions 
from the bunches of modes k and k' for which the phase 
9 ~ 7r, which corresponds to the interaction of packets 
of x and Sip waves. In contrast to the process of rescat- 
tering, the annihilation process XuXw — * ^0 k '"50 k + k '- k " 
and the inverse process will be suppressed because the 
corresponding time integrals have very rapidly oscillat- 

1— — cos vast 

mg exponents e m * . 

The analysis of Eq. ( |l33| ) shows the hard compo- 
nent 5<j) with k ~ fc* can be generated only during the 
very short time intervals At* ~ fc^ 1 around zeros of 
the inflaton field. The soft component with momenta 
k -C fc* is generated all the time. Soft particles produced 
at 1 0| > 0* have very small momenta in the range of 
< k < to. It makes sense to talk about such particles 
as free 0-particles removed from the coherently oscillating 
field 4>(t) only at time intervals r 3> toT 1 . An estimate 



of the soft component from (133) at the beginning of the 
process is (4> 2 } so ft — 9 2n \l rn %i whereas for the hard com- 
ponent one has (</> 2 )hard — (0 2 )soft/v / 9- Since 6(f) grows 
very fast, one h as to be careful with the range of validity 
of the solution (132). Indeed, Eq. (132) is only the first 
term in the iterative solution of Eq. (105). As soon as 6<j> 
grows , we have t o co nsider the iterative solutions of both 
Eqs. (104) and ( |105| ). We have to take into account the 
corrections to Xk due to the X and tp coupling on the 
r.h.s. of Eq. (104) as well as the second bilinear term on 



of Eq. 

r.h.s. ofEq. (105). Due to the exponential growth in the 



numb er o f particles, these corrections to the simple solu- 
tion (132) very quickly become important, which makes 
further investigation rather complicated. 

One should note that in addition to rescattering, there 
may exist other mechanisms of 0-particle production. 
For example, let us consider fluctuations 5<j) with effec- 
tive mass squared g 2 (x 2 )- As we already emphasized, 
this term is time-dependent. First of all, it experiences 
quasiperiodic high-frequency modulation, which, as we 
already noted in Sec. 
tional source of 1 
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may serve as an addi- 
-particles. In addition, the term g 2 (x 2 ) 
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oscillates with period During each oscillation it 

changes from its minimal value g&n x (for \4>{t)\ = $) 
to a much greater value ~ Sg^n^ 1 ^ 4 (for \4>(t)\ = 0*). 
This leads to a significant periodic change in the proper- 
ties of </>-particles, which is especially pronounced when 
1 4>(t) | A preliminary investigation of this issue 

indicates the possibility of a parametric resonance with 
^-particle production. 

Our main purpose here was not to give the final analy- 
sis of this issue but rather to outline different approaches 
to the problem of rescattering and ^-particle production, 
which should provide a proper framework for future in- 
vestigation. 



C. Rescattering and the end of preheating 



Can rescattering kill the resonance? In Sec. X A we 



found that rescattering can be rather efficient at the sec- 
ond stage of preheating. What can we say about the 
influence of rescattering on the development of paramet- 
ric resonance? 

The simplest idea would be to estimate the effec- 
tive mass of the %-particles induced by the fluctuations 
(4> 2 ): Am 2 , ~ g 2 (<fi 2 ). However, this would not be 
quite correct. Indeed, the whole process of x-particle 
production occurs in the interval \<f>\ < <f>„ during the 
time t* ~ (#m0$)~ 1/2 = k^ 1 , see Eq. @). If oscilla- 
tions of the modes 8<p occur during a longer time, then 
from the point of view of the creation of ^-particles they 
cannot be distinguished from the oscillations of the field 
(j){t) , and therefore they do not harm the development of 
stochastic resonance. We called such modes "soft," and 
the modes with k > fc*/4 "hard." 

Fluctuations of the scalar field <f> can be harmful to the 
development of the resonance if they can considerably 
alter the motion of the field (j> in the interval \<p\ < 
The only fluctuations which can change the direction of 



their motion during the short time t* 



are the 



modes with k > 2irk :t ^> fc*. This effect does not seem 
to be very important. At the time when the homoge- 
neous mode <p(i) enters the interval \<p\ > 0*, it has a 
kinetic energy <fi 2 /2 ~ to 2 $ 2 /2. In order to alter the mo- 
tion of the field <f> the "hard" fluctuations S(f> should (oc- 
casionally) have comparable (and opposite) speed, and 
therefore they should have a kinetic energy comparable 
to m 2 <I> 2 /2. Thus, the resonance disappears only after 
the kinetic energy of (/>-particles with momenta k 3> fc* 
becomes comparable to the total energy of the oscillat- 
ing field 4>{t). This could happen only at the very end of 
preheating. 

However, there is another mechanism which may harm 
the resonance. Each mode Xk "probes" space on a length 
scale Al ~ 2irk~ 1 . If the field 5(f) is homogeneous on 
this scale, it acts as a homogeneous background for the 
mode Xk- On the other hand, if S(j) is inhomogeneous 



on this scale, then the field Xk has an integrated inter- 
action with all inhomogeneities of the field 8<j> on the 
scale Al ~ 2irk~ 1 1 i.e. it interacts with the contribution 
to (<fi 2 ) from the modes with momenta greater than k. 
This corresponds to the appearance of an "effective mass 
squared" Am 2 , ~ g 2 (<l> 2 ), but only the modes with mo- 
menta greater than k should be taken into account in this 
calculation. Thus, from the point of view of the develop- 
ment of parametric resonance, one can introduce a new 
notion of an effective mass squared Am x (k) ~ g 2 ((j) 2 }k, 
where the index k means that we take into account only 
the modes with momenta greater than k. 

If the effective mass squared Am 2 (A:) becomes greater 
than k 2 , the equation of motion for such modes Xk 
changes considerably. This effect kills the resonance for 
the mode Xk if Am 2 (k) becomes greater than the width 
of the resonance. The resonance for the leading modes 
with k ~ fc*/4 ends when Am 2 (fc*) ~ g 2 (</> 2 )hard becomes 
greater than fc*/4. 

The difference between the total value of (cf> 2 ) and 
(<^ 2 )hard = (4> 2 }k,/4 can be quite significant. The num- 
ber of (^-particles produced in each scattering is equal to 
the number of x _ P ar ticles, each ^-particle taking away 
some portion of the momentum k of the corresponding 
X-particle. If this portion is small, 5<p fluctuations corre- 
sponding to these particles have momenta much smaller 
than fc*/4. Therefore, they do not give any contribution 
to the effective mass Am 2 (k <~ fc*/4), so they do not hurt 
the resonance at such momenta. If in the first collision 
a x - P ar ticle with momentum k ~ fc*/4 gives a signifi- 
cant portion of its energy to a ^-particle, then it loses its 
energy, and in subsequent collisions it will produce only 
harmless 6<p fluctuations with k <C fo*/4. 

Thus, one may argue that if rescattering is efficient, 
the number of "hard" ^-particles produced by x-particlcs 
should be similar to the initial number of x-particles with 
momenta ~ fc*/4, i.e. n^ ard < n x , whereas the total 
number of ^-particles produced by rescattering may be 
much greater. At the second stage of reheating, when 
g 2 (x 2 ) ^ m \i one can use an estimate 



1 

2^ 



k 2 dk 



x /k 2 +g 2 ( X 2 ) 



(134) 



If 



is dominated by soft fluctuations with k 2 <§; 



g 2 (x 2 )> then at the second stage of the resonance one 
should expect a strong anticorrclation between oscilla- 
tions of {x 2 } and (S(j) 2 }. This prediction is in agreement 
with the numerical results of |3Cf |. 

Now let us concentrate on the "hard" fluctuations with 
typical momenta ~ fc*/4. They can hamper the reso- 
nance if they make the field x massive, with an induced 



effective mass squared Am 



g 



/ hard 



comparable 



to the square of the typical momentum of x-P ar ticl es 
k ~ fc*/4: 



5 2 (<5</> 2 )hard >gm^ /16 . 



(135) 
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Suppose that a fraction 7 of all energy —3— is trans- 
ferred to the kinetic energy 4j-((5</> 2 )hard of "hard" fluctu- 
ations, 



gm<& 
32 



(^ 2 )ha 



This gives 



(136) 



(137) 



Comparison of Eqs. (137) and (135) shows that 7 > 
1/256, i.e. the resonance may slow down and eventually 
terminate only when the oscillating field <\> transfers at 
least ~ 1/256 of its energy to the "hard" fluctuations 
<\>. The total energy of all ^-particles will be somewhat 
greater than that. These particles get their kinetic en- 
ergy from the kinetic energy of x-particles ~ ^-|f— (x 2 )j 
so one may expect that the resonance terminates only af- 



ter 



32 



(x 2 ) becomes greater than 



256 



This can 



only occur close to the end of preheating. Let \J (x 2 ),. and 
<I> r be the values of x-Auctuations and amplitude of the 
background field at the moment t r when the parametric 
resonance is terminated by rescattering. Taking into ac- 
count that at the second stage of preheating m| ~ g 2 (x 2 } 
one finds that at the end of preheating 



VV)r>#r/16 



(138) 



Note also that \J (x 2 )r ^ 3>n because this would cor- 
respond to the result which we obtained in Sec 



IX B 



neglecting rescattering. In our subsequent calculations 
we will use the estimate \J (x 2 )r ~ 10 _1 $ r . This value is 
somewh at sm aller than \J (x 2 ) 2 — $2 which we obtained 
in Sec. IX B neglecting rescattering. However, the dif- 



ference between these two values is in fact not very large 
because $ r > $2- 

We are going to find \J (x 2 ) r and 3v, which should re- 
place our previous estimates for \J (x 2 ) 2 and $2 at the 
end of the second stage neglecting rescattering. Again 
we will use energy conservation. At the end of the first 
stage the energy density was equal to the potential en- 
ergy density to 2 $ 2 /2 of the inflaton field plus the energy 
of its interaction g^\n x ~ m 2 $ 2 , where m is the bare 
inflaton mass. At the end of the resonance (at the sec- 
ond stage), with an account taken of rescattering, the 
kinetic energy of the x-particles remains small, so the 
whole energy ~ 3m 2 i <I> 2 /2 transforms to the energy den- 
sity of interaction between x _ particles and the field cf>, 
Px = g$> r n x ~ g 2 (x 2 )r^ 2 r ~ 10"V$r- Note that p x 
includes the energy of the oscillating scalar field 4>{t). 
Energy conservation implies that <I> r ~ 3.5\/m&i/ g ~ 

— 1/4 

2.5$x9i • However, $ r obviously cannot be greater 
than $1. This means that r escattering can terminate the 
resonance either if \J (x 2 )r ^S> 10~ 1( l>r, in which case we 



essentially recover the previous results of Sec. IX B, or if 



Qi ^ 10 2 . In the last case one has \J (x 2 ) 
which yields 



0.35, 



/ mj>i 
9 
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(139) 



This estimate should replace Eq. (127) derived without 
account taken of rescattering. In particular, for g = 10~ 2 , 
which corresponds to q = 10 6 , and qi ~ 550, we get 
\J (x 2 )r ~ 2.5 x 10 15 GeV. To compare this result to the 
result of |j(J one should note that the definition of qo in 
|]30f differs slightly from ours, so it is better to compare 
our results for a given g rather than for a given qQ. In 
particular, one should compare their results for go = 10 6 
with our results for g = 10~ 2 : x/(x 2 )r ~ 3 x 10 15 GeV. 
This result agrees, to within a factor of 2, with the results 
of the lattice simulation of |30| . 

One should not overemphasize the significance of this 
agreement. The theory of the last stages of preheat- 
ing is extremely complicated, and there are many points 
in which our rough estimates could be improved. One 
should remember also that we are discussing stochastic 
resonance, which is extremely sensitive to even minor 
changes of parameters, see the table in Sec. VI [13 . From 



this perspective it is even somewhat surprising that one 
can describe many features of this process by analytical 
methods with rather good accuracy. 

Strictly speaking, the condition which we derived does 
not imply that the resonance is completely terminated. 
The leading modes, which have been amplified from the 
very beginning, stop growing when the effective mass of 
the field x becomes greater than k ~ fc*/4. However, 
the sub-leading modes still continue their growth until 
the effective mass becomes greater than k*/2. This pro- 
cess is very inefficient, but (x 2 ) continues slowly growing 
for a while. Moreover, (x 2 ) may grow a little even when 
the resonance is completely terminated and new particles 
are no longer produced. Indeed, due to the decay of the 
field 4>(t), the effective mass of the x particles becomes 
smaller, and therefore (x 2 ) may become greater even if n x 
remains constant. These effects are not very significant, 
but they make it difficult to clearly recognize the end of 
parametric resonance by looking at the behavior of (x 2 )- 
That is why throughout this paper, alongside the disper- 
sion of the fluctuations which is studied in most papers 
on preheating, we use the number density of particles to 
investigate the resonance. 

An estimate of the density of x-particles at the end of 
the resonance can be obtained by multiplying (x 2 )r by 
g^ r - 16#^(x 2 )r- It is given by 



0.4.9- 1 / 2 ( pmM^n- 1 - L2m 
\ 9 



Mr, 



.3/2 



(140) 



It is useful to compare this number with the number 
of ^-particles in the oscillating field <p(t) which sur- 
vive the rescattering. To distinguish the particles <f> in the 
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oscillating field and the free 0-particles created by rescat- 
tering, we will denote the number of particles in the clas- 
sical field as . At the end of the resonance it is given 

by m0$^/2, where is the effective mass gW (x 2 )r — 
0.1g$ r . Meanwhile n x ~ g$r(x 2 )r ~ I0~ 2 g$~. There- 
fore, 



(141) 



Eq. (141) says that at the end of the resonance x~ 
particles need to rescatter only 10 times to destroy the 
coherent oscillations of the classical field, i.e. to decom- 
pose it into separate </>-particles. Therefore one may ex- 
pect that at the end of the resonance or very soon after it 
X-particles may destroy the classical field <fi(t) completely, 
in agreement with ^TJ. This means that the final stage 
of decay of the homogeneously oscillating classical scalar 
field in our model is determined not by resonance but by 
rescattering. 

The decay of the classical scalar field <j>(t) is not the 
end of the story, but rather the beginning of a new stage 
of reheating. As we pointed out in |EJ, it does not make 
much sense to calculate the reheating temperature at this 
stage of the process. Indeed, from the point of view of 
the energy stored in the ^-particles, it is not very im- 
portant whether it is in the form of cf> fluctuations or in 
the form of a coherently oscillating field <fi. According 
to our estimates, the kinetic energy of x-particles may 
constitute only about 10~ 2 of the total energy at the end 
of parametric resonance. This estimate may be too pes- 
simistic, but even if the true energy is much higher, the 
main fraction of energy after the end of the resonance re- 
mains stored in the energy of ^-particles, and the energy 
of their interaction with x-particles. The total energy of 
%-fluctuations at large t decreases as a -4 , whereas the 
energy of (^-fluctuations as well as the energy of the os- 
cillating field 4>(t) at large t decreases as a~ 3 . Even if 
the total energy of the oscillating field cf>(t) and of <fi- 
particles were very small after preheating, eventually it 
woul d ag ain dominate the energy density of the universe. 



Eq. (141 



gives us additional information: the number of 
^-particles after preheating is at least ten times greater 
than the number of x-particles. If these particles do not 
decay, they will always dominate the energy density of 
the universe, which is unacceptable. Therefore when pre- 
heating ends one should apply the elementary (perturba- 
tive) theory of reheating [^|J6) to describe the decay of the 
remnants of the classical oscillating field tp(t) as well as 
the decay of the large amount of (^-particles created by 
rescattering. We will return to the theory of this process 
in a subsequent publication [pi. 



XI. PRODUCTION OF SUPERHEAVY 
PARTICLES DURING PREHEATING 

One of the most interesting effects which may become 
possible during preheating is the copious production of 



particles which have a mass greater than the inflaton 
mass to. This question is especially interesting in the 
context of the theory of GUT baryogenesis, which may 
occur in a rather unusual way if superheavy particles with 
mass M a few times heavier than m can be produced pi) . 
Such processes are impossible in perturbation theory and 
in the theory of narrow parametric resonance. However, 
we are going to show that superheavy %-particles with 
mass M 3> to can be produced in the regime of a broad 
parametric resonance. 

In order to study this regime let us return to Sec. 
where we made a simple derivation of the width of the 
resonance band, see Eq. (|32j). The only modification 
which should be made to this equation in the case where 
the field x has a ^-independent mass m x (0) = M is to 
add it to k 2 on the l.h.s. of the equation: 



M 2 < ( 3 2 0to $) 2/3 



2 J.2 

9 9 



(142) 



As before, the maximal range of momenta for which par- 
ticle production occurs corresponds to <f>(t) = 0*, where 

( t > * ~ h"\J m \ i ■ ^he maximal value of momentum for 
particles produced at that epoch can be estimated by 



+ M 2 



2 



The resonance becomes efficient for 
gm^ > AM 2 . (143) 



Thus, the inflaton oscillations may lead to a copious pro- 
duction of superheavy particles with M S> to if the am- 
plitude of the field $ is large enough, 3$ > 4M 2 /m. 

However, in an expanding universe $ and to^ are time- 
dependent. One should not only have a very large field 
at the very beginning of the process; one should continue 
to have gm<& > 4M 2 until the end of preheating. 

During the second stage of preheating both and 
$ change very rapidly, but their product remains almost 
constant because the energy density of the field <f), which 
is proportional to to 2 $ 2 /2, practically does not change 
until the very end of preheating. Therefore it is sufficient 
to check that gm& > 4M 2 at the end of the first stage of 
preheating. One can represent this criterion in a simple 
form: 



( gjiM p ln _! 10 
V 3to 



12 , 



1/2 



(144) 



For example, one may take M — 2m and g w 0.007, 
which corresponds to go = 10 6 in the normalization of 
Rcf. p0|. In this our condition (144) is satisfied, and 



an investigation with an account taken of rescattering 
shows a relatively insignificant suppression of (x 2 ), ap- 
proximately by a factor of 3. Our investigation suggests 
that for g ^> 10~ 2 this process should not be suppressed 
at all. Eq. (144) shows that for sufficiently large g one 
can produce superheavy particles with M 3> m. For ex- 
ample, production of ^-particles with M = 10to is possi- 
ble for g > 0.065. 
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In fact, suppression of superheavy particle produc- 
tion may be even less significant. Indeed, the reso- 
nance becomes strongly suppressed if it occurs only for 



As a result, the condition for 



the efficient preheating (143) can be slightly relaxed: 
gm^ > 2M 2 . This small modification implies that 
heavy particle production is not strongly suppressed for 

Af < m^ 4 



( 2 gfl M p , - 

^ 3m 111 9 5M P 
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For M = 10m 



this leads to a rather mild condition g > 0.036. 

We conclude that at least in our simple model, the 
production of superheavy particles is possible. However, 
with an increase of g the total n umb er of produced par- 
ticles becomes smaller, see Eq. (14C). It would be most 
interesting to investigate this issue in realistic models of 
elementary particles and to apply the results to the the- 
ory of baryogenesis. 



XII. DISCUSSION 

In this paper we discussed the theory of preheating for 
the simple model of a massive inflaton field <j> interacting 
with another scalar field \- As we have seen, the theory 
of preheating is very complicated even in such a simple 
model. Our main purpose was not to answer all questions 
related to the theory of preheating, but to develop an 
adequate framework in which these questions should be 
investigated. 

In the beginning particle production occurs in the 
regime of a broad parametric resonance, which gradually 
becomes narrow and then terminates. If the resonance is 
narrow from the very beginning, or even if it is not broad 
enough, it remains inefficient. We have found that broad 
resonance in an expanding universe is actually a stochas- 
tic process. The theory of this process, which can be 
called stochastic resonance, or stochastic amplification, is 
dramatically different from the theory of parametric res- 
onance in Minkowski space. Therefore one cannot simply 
apply the standard methods of investigation of paramet- 
ric resonance in Minkowski space; it was necessary to 
develop new analytical methods for the investigation of 
stochastic resonance in an expanding universe. We have 
found the typical width of the resonance ~ fc*/2 and the 
typical rate of the exponential growth of the number of 
produced particles in this regime. An important feature 
of our formalism of investigation of the broad resonance 
regime is its robustness with respect to modification of 
the form of the effective potential. Our methods should 
apply not only to theories with the potential m 2 <p 2 /2, 
but to any potential V(<p) when the resonance is broad. 

One should note, that the main reason why broad res- 
onance has a stochastic nature is the expansion of the 
universe. In the conformally invariant theories such as 

the theory j(f) 4 + ^~<^> 2 x 2 with 3 2 ^> A the resonance is 
broad but not stochastic because expansion of the uni- 
verse does not interfere with its development p9[. In 



realistic theories where the inflaton field <fi has mass m 
the conformal invariance is broken and one could expect 
that the broad resonance becomes stochastic as soon as 
the amplitude of the oscillations of the field <fi becomes 
smaller than m/vA. Indeed, for $ < mj\f\ the reso- 
nance is described by the model of a massive inflaton field 
considered in this paper. A more detailed investigation of 
this question shows that in models with g 2 3> A the res- 
onance becomes stochastic even earlier, at $ < -4= " 

VA oAJWp 

In our investigation of preheating we took into account 
the interaction of the oscillating inflaton field <j> with the 
particles produced during preheating. We have found, 
in particular, that the correction to the effective mass 
squared of the oscillating field tj> is proportional to , 
and the equation of motion of the field 4> acquires a term 
~ 5 n x]0|- This term experiences quasiperiodic oscilla- 
tions with a very high frequency ~ 2g$, which do not 
much affect the motion of the field <f)(t) but may serve as 
an additional source of </>-particles. 

We have found that if the coupling constant g 2 in the 

interaction term 4r</> 2 x 2 is small (j<3x 10~ 4 ), the res- 
onance terminates at the stage when the backreaction of 
produced particles is unimportant. For larger values of 
g 2 the resonance terminates due to a combined effect of 
the growth of the effective mass of the inflaton field and 
rescattering, which in its turn increases the effective mass 
of x-particles, making them heavy and hard to produce. 
We made an estimate of the number of x _ P ar ticles pro- 
duced during preheating and their quantum fluctuations 
(x 2 ) with all backreaction effects taken into account. 

Traditionally, the only purpose of the theory of reheat- 
ing was to obtain the value of the reheating temperature. 
From this point of view the theory of preheating for the 
simple model which we studied in this paper does not 
change the situation. For g <C 3 x 10 -4 the total en- 
ergy density of produced particles is exponentially small. 
Similarly, it remains extremely small even for large g if 
X particles have mass M much greater than about 10m. 
In the case when M is small and g > 3 x 10~ 4 , the x~ 
particle production is very efficient. However, we have 
found that even in this case after preheating one has 
many more ^-particles than x-particles. If x-particles are 
massless, or if they can easily decay, their contribution 
to the energy density of the universe rapidly decreases. 
Therefore, after preheating the main contribution to the 
energy density of the universe is again given by the <p- 
particles. The only difference is that prior to preheating 
these particles constitute the oscillating classical inflaton 
field <f>(t), whereas after preheating they acquire various 
spatial momenta and become decoherent. Thus, as we 
already pointed out in (!]], it does not make much sense 
to calculate the reheating temperature immediately after 
preheating. One should study the subsequent decay of 
the ^-particles . The theory of this decay is described by 
the elementary theory of reheating |]|||| ■ So why should 
one study extremely complicated nonperturbative effects 
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which may happen at the stage of parametric resonance, 
if in the end they will not greatly change our old esti- 
mates of the reheating temperature? 

We believe that the investigation of nonperturbative 
effects in the very early universe is worth the trouble. In 
fact, the complex nature of this process makes it espe- 
cially interesting. Indeed, a few years ago the standard 
picture of the evolution of the universe included a re- 
markable stage of explosive expansion (inflation) in the 
vacuum-like state, which is responsible for its large-scale 
structure, and a rather dull stage of decay of the inflaton 
field, which is responsible for the matter content of the 
universe. The processes which could happen during the 
later stage were typically ignored. 

Now we see that the stage of reheating deserves a more 
detailed investigation. Explosive processes far away from 
thermal equilibrium could impact the further evolution 
of the universe. As we know, the appearance of baryon 
asymmetry requires the absence of thermal equilibrium, 
so it is only natural to investigate the possibility of baryo- 
genesis at the stage of reheating, see e.g. ]|, 5]l[ll[ . 

Particles produced by the resonance have energies 
which are determined by the properties of the resonance 
bands. Typically this energy is much smaller than the 
temperature which would appear if the particles were in- 
stantaneously thermalized. Meanwhile, the total number 
of particles produced by parametric resonance is much 
greater than the number of particles in thermal equi- 
librium with the same energy density. Fluctuations as- 
sociated with these particles can be anomalously large. 
For example, we have found that for c ertain values of 
coupling constants in our model \J (\ 2 ) may become of 
the order of 10 16 GeV, and ^JW) may become even 
greater than W (x 2 )- In models describing several inter- 
acting scalar fields such anomalously large fluctuations 
may lead to specific nonthermal phase transitions in the 
early universe on the scale of 10 16 GeV [p|JlC|l. As we 
pointed out in j9|, the investigation of such phase tran- 
sitions in the theory of a single self-interacting field <j) is 
rather involved because one needs to separate the effects 
related to the oscillations from the effects related to the 
fluctuations of the same field. Therefore an optimal way 
to study nonthermal phase transitions is to investigate 
the models where the fluctuations produced during pre- 
heating restore symmetry for the field which does not 
oscillate during the oscillations of the inflaton field, see 
Sect. [X. We will return to the discussion of this effect 



in a separate publication J43J . 

Unlike fluctuations in thermal equilibrium, the non- 
thermal fluctuations produced by a parametric resonance 
often exhibit a nongaussian nature. In particular, "fluc- 
tuations of fluctuations" can be very large. This means 
that in some regions of the universe one can find fluc- 
tuations at a level much greater than its average value. 
This effect may play an important role in the theory of 
topological defect production. Indeed, even if the aver- 
age level of fluctuations is smaller than the critical level 



which leads to monopole production, they may be pro- 
duced in the rare islands where the level of the fluctua- 
tions is anomalously high. Note that in order to avoid 
cosmological problems and burning of neutron stars by 
the monopole catalysis of baryon decay, the density of 
the primordial monopoles should be suppressed by 20 to 
30 orders of magnitude. It was easy to achieve such sup- 
pression for the usual thermal fluctuations which appear 
after reheating, but for the nonthermal fluctuations pro- 
duced by resonance the situation may be quite different. 

There is an additional reason which makes the investi- 
gation of preheating so interesting. The theory of particle 
production in the early universe was one of the most chal- 
lenging problems of theoretical cosmology in the early 
70's. However, powerful methods of investigation devel- 
oped at that time produced rather modest results: parti- 
cle creation could be efficient only near the cosmological 
singularity, at densities comparable with M^. This pro- 
cess could not considerably increase the total number of 
particles in the universe. 

Now we see that in the context of inflationary cosmol- 
ogy all particles populating our part of the universe have 
been created due to quantum effects soon after the end 
of inflation. The investigation of these effects sometimes 
requires the development of new theoretical methods in- 
volving quantum field theory, cosmology, the theory of 
parametric resonance, the theory of stochastic processes, 
and noncquilibrium quantum statistics. 

In a situation where nonperturbative effects play an 
important role, and the number of produced particles 
grows exponentially, one could expect that the only re- 
liable tool for the investigation of preheating would be 
numerical simulations. Fortunately, one can go very far 
by developing analytical methods. For sufficiently small 
values of the coupling constant (g < 3 x 10~ 4 ) these 
methods allow us to make a very detailed investigation 
of preheating. For higher values of the coupling constant 
one can describe preheating analytically during most of 
the process. At the last stage of preheating the situa- 
tion becomes too complicated, and numerical methods 
become most adequate. Even in these cases analytical 
methods allow us to obtain estimates of the same order 
of magnitude as the results of numerical calculations, and 
sometimes this agreement is even much better. Taking 
into account all of the uncertainties involved in the an- 
alytical investigation of stochastic resonance as well as 
in the computer simulations, this agreement looks rather 
encouraging. It remains a challenge to develop a com- 
plete analytical theory of preheating, and to apply it to 
realistic inflationary models with many interacting fields. 
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